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1. Introduction 

Our main purpose is to present some recent topics for branching processes and 
superprocesses, with special focus on random media. The choice of the topics 
of course reflects the author's interest in the field and is necessarily incomplete. 
There are several excellent expository books on superprocesses in general, for 
example by Dawson, Dynkin, Etheridge and Perkins and here we suffice with a 
subjective choice of some recent developments. 

An outline of these notes is as follows. After defining the basic notions in 
Section [21 we turn to the connection between spatial branching processes and 
nonlinear partial differential equations (Section[3|). In the following section (Sec- 
tion |4|) we present some results on local extinction and local exponential growth, 
which is also a good excuse to discuss a version of the 'spine technology', a rel- 
atively recent (and, in the author's opinion, quite beautiful) tool appearing in 
the literature. The problem of the Law of large Numbers for spatial branching 
processes and superprocesses is discussed in Section [5l At the end of the 'de- 
terministic environment part', we picked two further topics in superprocesses: 
the so-called compact support property and the 'polar decomposition' for su- 
perprocesses. This second one is somewhat related to the Law of large Numbers 
discussed in the previous section, because when dealing with LLN, one many 
times first decomposes a certain 'weighted' (or /i-transformed) superprocess into 
its total mass ('radial' part) and a Fleming- Viot-type process ('angular' part). 

In the second part (Sections [7l fTT|) we will discuss spatial branching processes 
in random media, first giving a short review on some related classical problems 
(without branching), namely on Brownian motion among Poissonian obstacles 
and random walks in random environment (Section[7|). In the next three sections 
we introduce branching processes and superprocesses among Poissonian obsta- 
cles: hard obstacles, 'mild' obstacles, and some generalizations (Sections [3 O 
and[Tni respectively.) Finally, in the last section we review some basic catalytic 
branching problems. 

References (including work cited as well as further bibliographic items) are 
given at the end of each section. These are not intended to be full bibliographies 
on the subjects, but they are sufficient to give the backgrounds for the particular 
theorems and problems discussed. 

Acknowledgement. This paper is based on notes written for a graduate course 
at the University of Bath in July 2007. During the course several people, includ- 
ing A. Cox, S. Harris, A. Kyprianou, P. Morters, J. Pardo Millan and A. Sakai, 
contributed with valuable remarks to the lectures. The author owes thanks to 
each one of them. 

Several helpful suggestions by an anonymous referee are also gratefully ac- 
knowledged. Finally, the author is grateful to Prof. Aldous for inviting these 
lecture notes for publication in Probability Surveys. 
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2. Branching processes and superprocesses — basic notions 
2.1. Branching diffusions 

Spatial branching processes have two main ingredients : the spatial motion and 
the branching mechanism. We start with describing the spatial motion. 

Let L be a second order elliptic (differential) operator on the Euclidean do- 
main D C M'^ of the form 

l,J — l 2—1 

where aij,bi, i,j = l...,d, are in the class C^'^{D), rj e (0,1] (i.e. their first 
order derivatives exist and locally Holder-continuous), and the symmetric matrix 
{aij{x)) is positive definite on D. 

Of course, L can also be written in the slightly different 'divergence form': 

The operator L then corresponds to a diffusion process (or diffusion) on D 
in the following sense. Take a sequence of increasing domains Z3„ t D with 
-Dn CC I?n+]Q and let td^^ denote the first exit time from the (open) set D„. 
Then there exists a unique family of probability measures P^^ x d D {x is the 
starting point) describing the law of a Markov process Y satisfying that 

1. P^Yn^x) = 1, 

2- fiYtAToJ- Io''''""{Lf){Ys)ds is a martingale for all / G C^{D) and aU 
n>l. 

We say that the generalized martingale problem on D has a unique solution and 
it is the law of the corresponding diffusion process or L-diffusion Y on D. 

Note that it is possible that the event lim„_>co < oo ('explosion') has 
positive probability. In fact lim„^oo td^ < oo means that the process reaches 
A, the 'cemetery state' in finite time, where A is identified with the Euclidean 
boundary of D in the case of a bounded domain, or the Euclidean boundary of D 
augmented by a point 'at infinity' in the case of an unbounded domain. In other 
words, the process actually lives on D ^ DU{A}, the one-point compactification 
of D and once it reaches A it stays there forever. 

The connection to linear partial differential equations is well known. For a 
bounded continuous function /, consider the parabolic equation: 

it — Lu in (0, oo) x D (1) 

with initial condition 

ii(0, •) = /(•) mD. 



the notation A GC B means that A is bounded and its closure is still in B. 
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This Cauchy problem (generalized heat equation) is then solved by u, where 
u{t,x) :— Exf{Yt) and Y is the diffusion corresponding to L on D. 

Furthermore, if Ttf is defined by Tt{f){x) := E^f{Yt), then the Markov 
property of Y yields that {Tt}t>o is a semigroup, that is Tt+s = Tt o Tg for 
t,s>0. Now, dl]) gives 

Tt+hf{x)~Ttf{x) rirrf\(\ 
hm = L[Ttf)[x), t > 

hlO t 

and formally we obtain {t = 0) that 

hm Lf, 
hio t ■' 

which can in fact be verified for a certain class of f s. Hence L is often referred 
to as the infinitesimal generator of Y . 

Diffusion processes behave in a very nice way from the point of view of their 
large time behavior, namely there are exactly two cases. Either 

Vx e -D, 7^ -B CC D open, Px{Yt £ B for arbitrarily large t's) = 1 

or 

yx e D,(l) ^ B dd D open, P^iYt ^ B for all t > T{B,uj)) = 1. 

In the first case we say that Y is recurrent and in the second we say that Y is 
transient. 

After describing the motion component (L-diffusion) we now turn to the 
branching component. Suppose that we start with a single individual and any 
individual has 0,1,2,... offspring with corresponding probabilities pQ,pi,p2, ■■■ 
and suppose that branching occurs at every time unit. If h is the generating 
function of the offspring distribution, 

h[z) :=po +Piz+p2Z^ + 

then clearly 

/i'(l)=Pi + 2p2 + 3j33 + --=:e, 

and e is the expected (or mean) offspring number. 

As is well known, there are three cases for the Galton- Watson process defined 
above: 

1. /i'(l) < 1 and the system dies out in finite time [suhcritical case) 

2. /i'(l) = 1 and the system dies out in finite time {critical case) 

3. /i'(l) > 1 and the system survives with positive probability {supercritical 
case) 

For example we are talking about strictly dyadic branching when P2 = 1 and in 
this case h{z) — and h! (\) — 2. 

In the supercritical case e„ , the expected offspring number at time n, satisfies 
tn = e" and when po > 0; = P(extinction) is the only root of h{z) = z in 
(0, 1). (The function h can be shown to be concave upward.) 
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Many times, instead of unit time branching one considers random branching 
times with exponential distribution. In this case the above classification of sub- 
critical, critical and supercritical branching is the same as before. Moreover, if 
the exponential rate is /? > 0, then the probability of not branching up to time 
t is obviously e~^* and a well known limit theorem due to Biggins states that if 
Zt is the population size at t then in the supercritical case 

3 lim e-'^^Zt = N >0, a.s., 

t — >oo 

on the survival set, as long as the offspring number X satisfies that XlogX has 
finite expectation. (There is a similar theorem for unit time branching too.) 

After discussing motion and branching separately, let us try to put the two 
ingredients together! To this end, let /3 > be smooth (more precisely. Holder 
continuous with some parameter rj £ (0, 1]) and not identically zero. The (L, /3) — 
branching dijfusioi^ X has the following (informal) definition. 

Definition 1 {{L,j3)— branching diffusion). Start with one initial particle lo- 
cated at some point x Cz D, performing an L-diffusion on D. Her exponen- 
tial clock depends on her position in the sense that the instantaneous rate of 
branching at y € D is /3{y). So for example the probability the particle has not 
branched by i > is exp(— /3(Kj) ds). When the clock rings she is replaced 
by two offspring (strictly dyadic branching) The two offspring, starting from 
the position of the parent particle perform independent L-diffusions and are 
equipped with their independent exponential branching clocks, where again the 
instantaneous rate of branching at y G D is P{y); etc. This way one obtains 
the stochastic process Z (the value of Zt, t > is a particle configuration, or 
discrete measure on D); we will denote by Zt{B) the number of particles at t 
in BQD. 

Since (3 is spatially dependent it is no longer clear what one means by sub- 
critical, critical or supercritical case. Also, the criterion for extinction becomes 
a nontrivial question. One has to distinguish between extinction and local ex- 
tinction. 

Definition 2. Z exhibits local extinction if for all B CC D, there exists an a.s. 
finite random time tb such that Zt{B) = for all t > tb- 

Intuitively, whether or not local extinction occurs should depend on the trade- 
off between two properties: 

1. how transient Y (corresponding to L) is, and 

2. how large f3 is. 

So, for example if Y is 'very transient' comparing to the size of f3, then the 
motion component 'wins' and Z exhibits local extinction. 

^Sometimes, when the domain is emphasized, one writes {L, /3, D)— branching diffusion. 
^Obviously one could consider more generally a random number of offspring (possibly zero) 
according to any given branching law. 
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The 'strength of transience' can be measured by a single number Xc{L) < 
called the generalized principal eigenvalue of L on D. We will give a different 
definition later but note here that for i? CC -D open and x G B, 

lim - \ogPa,{Ys e B, for aU s < t) = Xb < 0, 

t— >oo t 

where Xb does not depend on x £ B, and Xb„ 1 Xc{L) for a sequence of domains 
B„ t D. So Ac is loosely speaking the 'exponential rate of escape from compacts'. 

The main result concerning local extinction (we will discuss its more general 
form later) is that when /3 is a positive constant, local extinction occurs if and 
only if Ae(L) < -/3. 

We close this section with two definitions of superprocesses — as scaling 
limits of branching diffusionf|f| and also as measure valued Markov processes via 
their Laplace functionals. We start with the latter one. 

2.2. Superprocess via its Laplace functional 

Let D C K'' be a domain and let B{D) denote the Borel sets of D. We write 
A4 f{D) and A4c{D) for the class of finite measures resp. the class of finite mea- 
sures with compact support on B{D). For € Mf{D), denote := fJ.{D) 
and let C^{D) and C^{D) be the class of non-negative bounded continuous 
resp. non-negative bounded continuous functions D ^ M. having compact sup- 
port. Write C'^'^{D) for the usual Holder spaces of index rj £ (0, 1] including 
derivatives of order k, and set C"'(L') := C°'''{D). 
Let L be given (in divergence form) by 

L iy -aV + fe- V, (2) 

where aij,bi £ C^'^{D), i,j = l,...,d, for some rj £ (0,1], and the matrix 
a{x) :— (ai.j(x)) is symmetric, and positive definite for all x £ D. In addition, 
let a,(3 £ C^{D), and assume that a is positive, and (3 is bounded from above. 

Definition 3 ((i, /3, a; -D)-superdiffusion). With D,L,(3 and a as above, let 
{X,Pfj^ , fi £ Aif{D)) denote the (L, a; Z?)-superdiffusion, where /i denotes 
the starting measure Xq. That is, X is the unique A^y(Z?)-valued continuous 
(time- homogeneous) Markov process which satisfies, for any g £ C^{D) , 

Ef,exp{Xt,-g) = exp{fi,-u{-,t)), (3) 

where u is the minimal nonnegative solution to 

ut = Lu + (3u — av? on D X (0, oo), 

l™w(-,t) = g(-). 

As usual, {v, /) denotes the integral v{dx) f{x). 

*As already mentioned, branching diffusions are also measure valued processes: Zt{-) = 
X^i^l ''z»(')' where Nt is the population size at time t. 
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One usually refers to (3 as mass creation and a as the intensity parameter 
(or variance). Finally, the Markov property is in fact equivalent to the property 
that the shift defined by (|4]) defines a semigroup, which in turn follows from the 
minimality of the solution. 



2.3. The particle picture for the superprocess 



Previously we defined the (L, f3, a; D)-superprocess X analytically. In fact, X 
also arises as the short life time and high density diffusion limit of a branching 
particle system, which can be described as follows: in the n}^ approximation 
step each particle has mass 1/n and lives a random time which is exponential 
with mean 1/n. While a particle is alive, its motion is described by a diffusion 
process corresponding to the operator L (on D — D U {A}). At the end of its 
life, the particle located at x € D dies and is replaced by a random number of 
particles situated at the parent particle's final position. The law of the number 
of descendants is spatially varying such that the mean number of descendants 
is 1 + while the variance is assumed to be 2a{x). All these mechanisms are 
independent of each other. 

More precisely, for each positive integer n, consider Nn particles, each of mass 
i, starting at points x^"^ £ D,i = 1,2,...,7V„, and performing independent 
branching diffusion according to the motion process Y , with branching rate 
cn, c > 0, and branching distribution {Pfc"''(a;)}^Q, where 



e„(x) = fcp^"-' (x) = 1 

fc=0 

and 



vl{x)^Y.i^~lfpt\x)=m{x)+o{l) 



k=0 



as n — > oo, uniformly in x; m, 7 G C^{D), rj € (0, 1] and m{x) > 0. Let 



n -"^ — ' 

i=l 



Let Nn{t) denote the number of particles alive at time t and denote their posi- 
tions by {Xl'(t)}fj:l*\ Denote by MfiD) {Mf{D)) the space of finite measures 
on D (D). Define an Ai f{D)- valued process Xn by 



n 

1=1 



Denote by P^"'' the probability measure on D {[0, 00), Mf{D)) induced by Xn- 
Assume that m and 7 are bounded from above. Let w — lim„^oo IJ-n = f^ & 
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A4f{D). (The notation w — lim denotes the hmit in the weak topology.) Then 
there exists a P* e C{[0,oo),Mf{D)) such that P* = w - hm f/,"^ Define 

Pf, e C([0,oo),7W/(-D)) by Fp(-) = n D) and let X be the process corre- 
sponding to Pfj_. Then X is an (L, /?, a; I?)-superprocess, where L corresponds 
to Y on D, f3{x) — c'y{x) and a{x) — ^cm{x). 

Finally, in the particular case when /3 = 0, an alternative approximation 
yields the same superprocess: in the n^^ approximation step one considers critical 
branching difi^usions (the motion component corresponds to L and the branching 
is strictly binary, i.e. either zero or two offspring), but the branching rate is now 
2na{x). So a(-) in this case can also be thought of as the branching 'clock'. 
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3. Connection between spatial branching processes and nonlinear 
partial differential equations; the qualitative behavior of 
superprocesses 

3.1. Log-Laplace equations 

Let us recall the log-Laplace equation for superprocesses: for any g e C'^{D), 

E^exp(Xf ,-5) = exp {fi,~u{-,t)), 

where u is the minimal nonnegative solution to ([4]). This fact already shows 
that there exists a connection between measure valued processes and partial 
differential equations. 

Remark 1. A similar connection with the (L, /3)-branching diffusion Z (with 
probability P) is as follows. For any g S C'^{D) , u{x,t) = E^exp {Zt ,—g) 
solves 

ut = Lu + (3(v? ~ u) on D X (0, 00), 
limu(-,i) = e"9(-)^ 

< li < 1. 

Note that, unlike in the superprocess case, here the term simply comes from 
the fact that the branching is strictly dyadic and so the generating function is 
just h{z) = . Note also that in ([5]) one can always 'switch' from — u to 
V — as w 1 — u solves 

vt ^ Lv ^ (3{v ~ v^) on £)x (0,00), ' 
limi;(-,0 = 1 - e~»(-\ 

< w < 1. 

o 

3.2. Exit measure and Brownian snake 

Dynkin's exit measur^ yields yet another connection with semilinear partial 
differential equations. Informally, the exit measure from a domain Z? C M'' 
is obtained by 'freezing' mass of the superprocess when it first hits the boundary 
dD. One way of defining the exit measure is through the approximating particle 
picture: the exit measure can easily be defined for the nth level discrete particle 
system, and then one has to verify that as n ^ 00, those discrete exit measures 
converge to a limit. (See [5] for more on exit measures.) 

Although in this subsection we will talk about Brownian motion, in fact the 
method can be extended from Brownian motion (and equations involving A) to 
general diffusions (and equations with elliptic operators L). 




""We will revisit a version of it in Subsection 14. II 
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Dynkin showed that if D is smooth and 5 is a bounded function on 9D, then 
u{x) = — logMsx cxp—{X^,g),x e D solves the problem 

Au = on D, ~\ 

u\dD = 9- J 

In order to investigate solutions with boundary blow-up, choose gn = n and then 
let n t 00, quickly leading to the following representation^ 

Corollary 1. Let D cW^ be a bounded regular^ domain. The function u{x) = 
— logPsx{X^ = 0),x E D is the minimal nonnegative solution to the problem 

Au = on D, ~| 

u\aD = +00. J 

Before turning to the next part, let us quickly mention that Dynkin actually 
defined the exit measure not only from spatial but also from space-time domains 
(c.f. Subsection 14. ip . leading to probabilistic solutions to nonlinear parabolic 
boundary value problems. (We do not discuss this further here, the interested 
reader should consult [5].) 

Le Gall introduced a path-valued process called the Brownian snake and 
related it to the exit measure, thus providing yet another probabilistic repre- 
sentation for nonlinear PDE's (in terms of the snake). 

The Brownian snake combines the genealogical structure of random real trees 
with spatial motions governed by a general Markov process (usually taken as 
Brownian motion). Informally, each Brownian snake path corresponds to the 
spatial positions along the ancestral line of a vertex in the branching tree. The 
precise definition of the Brownian snake is motivated by the 'coding' of real 
trees. In view of applications to PDE, Le Gall characterized the exit measure 
as a measure, which, in a sense, is uniformly spread over the set of exit points 
of the Brownian snake paths from D. (See for the precise definition of the 
snake and the characterization of the exit measure.) 

When the underlying spatial motion is d-dimensional Brownian motion, the 
connection between the exit measure and the semilinear PDE Au = u^ is as 
follows. 

Recalling that the snake is a path valued process and that the exit measure is 
expressed in terms of the snake, Le Gall defines a a-finite measure on C(R+, W) 
denoted by Njf|, where W is the set of all finite paths (see [S] for the precise 
definition of N^:), and proves the following representation (notice the extra factor 
4 in the equation). 

"in |12l [S] there is a different approach to solutions with boundary blow-up without using 
the exit measure. 

'^The domain D is called regular if every x € dD satisfies that inf{t > | ^ D} = 
Px-a.s., where (5,P) is Brownian motion 

*It is dubbed the 'excursion measure away from x'. 
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Proposition 1. Let D be a bounded regular domain. Let g be a continuous 
nonnegative function on dD. Then u{x) = Nj;(l — cxp — {X^,g)),x E D solves 
the problem 



u\dD = 9- J 

Similarly as before, this leads to the following representation of the solution 
with boundary blow-up. 

Corollary 2. Let D be a bounded regular domain. Then the function u{x) = 
l^x{X^ ^ 0),x E D is the minimal nonnegative solution of the problem 



Due to groundbreaking work of Dynkin, Le Gall and Mselati (see [S] |9l [TT]), 
we now have a good understanding of the classification and probabilistic repre- 
sentation of solutions of Am = in a domain. 

3.3. Support properties 

It turns out that (jl]) can be used very efficiently in investigating the following 
properties of the superprocess. 

Definition 4. A path X. survives if Xt ^ 0, Vt > and becomes extinct 
otherwise. 

We will use the shorthand S := {Xt ^0, Vi > 0} for the event of survival. 
The basic question is to calculate Pp(5), and in particular to decide whether 
this probability is zero ('extinction'). 

Definition 5. The superprocess exhibits local extinction if for every set B CC 
D there exists a P^-a.s. finite random time tb such that Xt{B) =0, Vt > tb. 

Unlike for discrete particle systems, however, it is not clear whether 'i? be- 
comes eventually empty' is the same as \imt^ oo Xt{B) = 0. R. Pinsky [12] 
introduced the following notions. 

Definition 6. Assume that P^{S) > 0. 

• The support of X is recurrent if 





Pf,{Xt{B) > 0, for some 



t\ S) ^1, WB C D open. 



• When d>2, the support of X is transient if 



V^iXtiB) > 0, for some t | 5) < 1, VB CC D s.t. 



D \ B is connected and supp(/i) O B — (d. [B is the closure of B.] 
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• When d — 1, the support of X is transient if 

P^,{Xt{B) > 0, for some t \ S) < I, 

either for aU B CC D s.t. inf supp(/i) > supS or for aU B CC D s.t. sup 
supp(/i) < inf B. 

It turns out that in the above definitions there is no dependence on the 
starting measure, i.e. each property (extinction, survival, local extinction, re- 
currence and transience) are either true for all 7^ /i G Ale or not true for any 



3.4- Weighted superprocess and nonlinear h-transform 

A very important tool in the investigations of these support properties is a 
transformation that leaves the support invariant but changes the mass. If X is 
a superprocess and ft. > is smooth, then we will define below a weighted super- 
process X^ which has the same support but has different motion and branching 
mechanism. Let 

dX^ 

X^ := hX, i.e — f =h, t>0, 
dXt 

and notice that the log-Laplace equation with (j) — hip, fi — v/h yields 

where u is the minimal nonnegative solution to (jj]) with in place of g. 

It is sufficient to restrict the setting on compactly supported measures, v G 

Mc <^ fJ-e Mc and test functions, (f) G C+{D) <^ ^ e C+{D). 
Let A{u) := Lu + (iu — au^. Now, if v :— u/h, then 

where ^ 

Vt = yA{hv) 

ri 

and v{-,0) = V'(')- What this means is that X'^ is another superprocess! If we 
now introduce the notation A'^{-) :— j^A{h-) (analogously to the Doob's h- 
transform for linear operators) , then we obtain the log-Laplace equation for the 
weighted superprocess X'^ in the following form: for any tp G C^{D), 

E^exp{Xt,-g) =exp{i^,-v{-,t)), (6) 

where v is the minimal nonnegative solution to 

Vt = A!'(v) on D X {0,00), 



limM(-, t) = V-'(-)- 



(7) 
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That is, if X corresponds to the operator A then X'' corresponds to the operator 
A'^. Although X''- is just a weighted version of X, it corresponds to a different 
motion and a different branching mechanism, since, as a straightforward com- 
putation reveals, 

A'^iv) = ( L + a— ■ V + -^^ j^-^ j V - ahv. 

So the new motion corresponds to L+a^ • V and the new branching mechanism 
corresponds to the term i^d-Slhy _ ahv^. 

Just like the transformation X X^ preserves the support of X, the cor- 
responding transformation A — > A''' (nonlinear /i-transform) preserves some 
important properties of the operator. For example, the generalized principal 
eigenvalue of the linear part of the operator is invariant under the transforma- 
tion. 

3.5. Some results 

In the previous section we discussed the generalized principal eigenvalue of a 
diffusion operator and its connection to local extinction for branching diffusions 
(with constant branching rate). We now discuss it in more generality. 
Define the generalized principal eigenvalue as 

Ac = Xc{L + P) :== inf{A | 3m > s.t. (i -f- /3 - A)m = 0, } 

and note that for constant (3, one has \c{L + l3) — \c{L) + (3. Now a more general 
result concerning local extinction is as follows. 

Theorem 1. The {L, (3, a, D)-superprocess (or the {L, (3, D) -branching diffu- 
sion) exhibits local extinction if and only if Xc < 0. 

(The superprocess version of this theorem is due to R. Pinsky.) For super- 
processes, since extinction implies local extinction, one obtains that a sufficient 
condition for survival w.p.p. is that Xc{L) > 0. 

It is beyond the scope of these notes to discuss all the support properties, 
but as an illustration, we show how one can handle the extinction probability 
by using the log-Laplace equation. 

First, notice that for jj, £ A4c, 

PJXt = 0) = E„ lim e<-"'-^'> = lim E„e<-"'^*> = lim e<-"'"'(-'*)'''> , 

where is the minimal nonnegative solution to ([4]) with g replaced by n. It is 
clear probabilistically that u'"^ is monotone increasing in n and it can be shown 
that u := lim„^oo is finite and it is the minimal nonnegative solution to ([4]) 
with g replaced by oo. 
We have obtained 

P;.(Xt = 0) = e<-<'*)'^>. 
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and now letting t "f oo and defining w{-) ~ linit^oo i) (again, monotonicity 
is clear probabilistically), we obtain the probability of extinction: 

As expected, it can be shown that w solves the steady-state equation Lw + (3w — 
aw^ = 0. The question is now: which (nonnegative) solution of the steady-state 
equation is w7 This is a subtle issue as w is sometimes the maximal solution 
and sometimes it isn't. 

More precisely, it can be shown that the maximal solution Wmax always exists 
and 

P^(7^ CC D) = e^-"-— 

where TZ Us>o supp(-'^s) is the range of the process. 
Since w < uTmax, Pp(7^ CC D) < P^(5^). 

Suppose now that the superprocess possesses the compact support property, 
that is, if fi is compactly supported then so is 1Jq<-^<( supp(Xs) for all t > 
(this property too is independent of the starting measure). Then, obviously 
P^(7^ CC D) > P^(5^), that is P^(7^ CC D) = P^(5^) and w = w^ax- 

This observation leads to a nice way of proving explosion (no compact support 
property) for the superprocess. Namely, all one has to show is that w is not the 
maximal nonnegative solution to the steady-state equation. For example, when 
we know that the process becomes extinct (and thus w — 0), the existence of a 
positive solution yields explosion. 

Based on these ideas one can come up with examples where the underlying 
process is conservative (i.e. no explosion occurs), but nevertheless the super- 
process explodes. This can be shown 7 to be true even in the case of a time 
changed one dimensional Brownian motion with /3 = 0, a = 1 (i.e. in the case of 
a (/3(a;)AM, 0, 1; ]R)-superprocess), if the time change p > is sufficiently large. 
(Note that the time change does not effect recurrence and so the time changed 
Brownian motion is still conservative.) 

Of course one may consider the {p{x)Au, 0, 1; ]R'*)-superprocess in any dimen- 
sion d > I, and wonder, how large p has to be in order to loose the compact 
support property. Somewhat surprisingly, it turns out ^ that while the thresh- 
old for p is quadratic for d > 2, it is cubic for d = 1. 
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4. Local extinction versus local exponential growth; the 'spine' 

In this section we present a 'spine' construction and its application in the study 
of local extinction versus local exponential growth of spatial branching processes. 
Although it can be applied to superprocesses too, we will only discuss, the more 
intuitivqj branching diffusion case. 

Let X denote the {L, /3)-branching diffusion in the section. The finite (dis- 
crete) starting measure will be denoted by /i. 



4.1. The spine decomposition 



A 'natural martingale'. Let B CC D be a nonempty domain. Suppose that 
supp (fi) C B. Let J^t denote the natural filtration of X up to time t > and 
let X*'^ denote Dynkin's exit measure from B x [0,t). Loosely speaking, X*'^ 
is obtained by 'freezing' particles on dB at the time they hit it and so X*'^ is 
a discrete measure on {B x {t}) U {dB x (0,i)). 

Let A denote the generalized principal eigenvalue^"! of L + /3 on B. It is a 
standard fact that there exists a unique nonnegative function (f> which solves 
{L + f] — X)(j> = on B with Dirichlet boundary condition. Define for each fixed 
^> 0, 0* : B x [0,t] [0,00) such that </>*(•,«) = </>(•) for each u G [0,t]. (Here 
B denotes the closure of B). Wc claim that 

Mf:^e-^'^-^^/^^:t>0 

is a (mean one) P^-martingale. To see this note that on account of the so- 
called Dynkin's Markov property of exit measures and the Dirichlet boundary 
condition for 0, 



t+s -^t+s,B 



Tt I = e-^'Ex, I e-^^ 



,X' 



for s, i > 0. We have to show that the righthand side equals Mf a.s. On account 
of the branching property it is enough to show that Eg^ (e~^* ((/)*, X*'^)) — (j){x) 
for all t > and x E B. To show this latter property, note that from the 
stochastic representation for backward equations: v{x, 0) :— Eg^ (e"'*'* ((/)*, X*'-^)) 
where v solves —v = {L + (3 — X)v on B x (0, t) with v — cj)^ on 9*'^ and 9*''^ 
is the union of dB x {0,t] and B x {t}. The boundary condition follows from 
the regularity properties of the underlying diffusion. Therefore, by parabolic 
uniqueness, we obtain that f(-,0) = </>(•); t>0. 

Our goal is to define a change of measure with this 'natural' martingale. (For 
the reader unfamiliar with this notion: if the new measure Q for a continuous 



®The author's wife is not a mathematician. After hstening to phone conversations with A. 
Kyprianou she asked 'who this Martin Gale was' and 'what was wrong with his spine'. 



^''which is just the usual principal Dirichlet eigenvalue. 
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time stochastic process with fihration {Gt} is defined by 

= Mt, (8) 



St 

then it is easy to show that M must be a mean one P^^-martingale. Conversely, 
every mean one martingale defines a proper law for the stochastic process with 
the definition under ([8]).) Before we do that, let us first define two changes of 
measures: the first is for the motion component and the second is for a Poisson 
point process. 

As before B will always denote a nonempty open set compactly embedded in 
D with a smooth boundary. 

Girsanov change of measure. Let X = Xc {L + /3, B) . Just like before, </) 
denotes the eigenfunction satisfying 

(L + /3 - A) = in B with = on dB. 

Let = 'mi{t > : Yt ^ B} and assume that the diffusion (Y.Fx) is adapted 
to some filtration {Qt : t > 0} . Then under the change of measure 



Qt 



exp<^-y X-P{Ys)ds 



the process (F, P^) corresponds to the /i-transformed [h — 0) generator (L + 
/3 - A)-^ = L + a(/.-iV0 • V. 

It can be shown that the process {Y, PJ^) is ergodic on B (i.e. it is positive 
recurrent). 

Change of measure for Poisson point processes. Given a nonnegative 
bounded continuous function g{t), < > 0, the Poisson point process rj is defined 
as follows. Let {n = nt;t > 0} be a Poisson process with instantaneous rate 
g{t), and let rj — {{ci : i = 1, rit} : t > 0} be the corresponding Poisson point 
process on [0, oo) {ai is the ith 'arrival'). Let denote the law of rj. 

If rj is adapted to {Qt : t>0}, then 



dha 



Qt 



2"'exp|-^*g(s)ds| 



Combining the change of measure for the motion component and for the 
Poisson point process leads to a change of measure for the spatial branching 
process X . Our 'spine' result here is as follows. 

Theorem 2. Suppose that fj, is a finite measure with supp/i C B. For branching 
particle process we can thus take fi = ^^5xi where {xi\ is a finite set of (not 
necessarily distinct) points in B. Define by the martingale change of measure 
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Define 

that is, we randomize the starting point of (YjFt) according to the probability 
distribution (j)fi/ {(p, fj) . Note in particular that when ^ — 6x, — P^- 

Suppose that g G C^^ (D) and Ug is the minimal positive solution to u = 
Lu + /3u'^ — (3u on D X (0, oo) with lim^io t) = 9i')- Then 

E,[e-<^'^^^) = 

Y: I EtL^^(^^ (e-^(^') f[ug{Y^^,t-a,))ll Ug{x„t) \ (9) 

i [ k=i j^i J 

The decomposition suggest that (X, P^) has the same law as a process con- 
structed in the foUowing way. From the configuration = pick a point 
x' G {xi\ with probabihty <j)[x') / {(j)^ /j). From the remaining points, independent 
(L, /3; D)-branching processes are initiated. From the chosen point, a (Y,P^,)- 
diffusion is initiated along which (L, /3; _D)-branching processes immigrate at 
space-time points {{Ya-. , ai) : i > 1} where n = {{ai : i — 1, ut} : i > 0} is a 
Poisson process with law L^''*^^^ 

This decomposition relates to a spectrum of similar results that exist in 
the literature for both superprocesses and branching processes by S. Evans, 
A. Etheridge, R. Williams, S. Roelly, A. Rouault and others 0- Theorem E] of- 
fers decompositions with the particular feature that the spine is represented by 
a diffusion conditioned to stay in the compactly embedded domain B. 

Now, obviously the new measure is absolutely continuous with respect to the 
old one up to time t. However this does not mean that this is true 'up to time 
infinity'. Nevertheless, we have the following result. 

Lemma 1. Suppose that supp/i C B and A = Xc{L + f3,B) > 0. Then Mf 
converges to its almost sure limit in (Pfi) , and furthermore P^^ ^ P^,. 



4-2. Application to local extinction and exponential growth 

After all these preparations, let us see how one can attack the problem of local 
extinction and (local) exponential growth for the branching diffusion. 

Theorem 3 (Local extinction vs. local exponential growth). Let Q ^ fi be a 

measure with supp/i dd D. 

(i) Under P^ the process Z exhibits local extinction if and only if there exists 
a function h > satisfying {L -|- /3)/i — Q on D, that is, if and only if 
Ac < 0. In particular, the local extinction property does not depend on the 
starting measure fi. 

^^The spine decomposition presented in this notes was proved by A. Kyprianou and the 
author and to a large extent it was inspired by a result of R. Lyons, R. Pemantle and Y. Peres. 
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(ii) When Ac > 0, for any A < Ac and any open % ^ B CC D, 

\ limsupe"^*Zt(B) = 00) >Q andP^i limsupe"^^*Zt(B) < 00 ) = 1. 

Remark 2 (Total mass). In Theorem|3]we are concerned with the /ocaZ behavior 
of the population size. When considering the total mass process \\Zt\\ := Zt{D), 
the growth rate may actually exceed Ac- Indeed, take for example the (L, /3; D)- 
branching diffusion with a conservative diffusion corresponding to L on Z? and 
with Ao := Xc{L,D) < 0, and let /3 > be constant. Then Xc{L + /3,D) = 
P + Xq < (3, but - since the branching rate is spatially constant and since there 
is no 'loss of mass' by conservativeness - a classical theorem on Yule's processes 
tells us that e~^*\\Zt\\ tends to a nontrivial random variable as t ^ 00, that is, 
that the growth rate of the total mass is (3 > Xc- o 

Before the proofs we present a proposition, which is not hard to prove by an 
application of the Borel-Cantelli lemma. 

Proposition 2. For any nonempty open set B <Z D and finite fj,, 

P^ (limsupZt(B) e {0,c»} ) - 1. 
V tToo / 



4.3. Proof of Theorem\3\ (i) 

(a) Assume that Ac < 0. Then there exists an ft, > solving (L + (3)h — 0. We 
claim that {(/i, Zt) : i > 0} is a positive P^-supermartingale for all supp(/.t) CC 
D. Indeed, it is not hard to show that Es^{h, Zt) < h{x) for t > and x G D. 
From a standard application of the Markov property and the branching property, 
it then follows that {{h, Zt) : t > 0} is a P^-supermartingale. 

Now in the possession of this supermartingale it follows for Borel B CC D, 
that 

Hm sup Zt (B) < C hm sup(/i, Zt) < 00 (10) 

tloo tlco 

P^-almost surely where C is a constant. When B is open, by Proposition [2] it 
follows that Hmfjoo Zt [B] = P^ — a.s. Since every compactly embedded Borel 
can be fattened up to an open B CC D, local extinction follows by comparison. 

(b) Assume now that Ac > 0. Since it is assumed supp/j, CC D, we can 
choose a large enough B for which supp /i C P and A = Xc{L + (3,B) > 0. 
Choose ^ g £ C+ so that g < Is; obviously, it suffices to prove that 
Pfi (limsupt|oc (g, Zt) > 0) > 0. Let m/ and P^ be as in Theorem [2] By Lemma 
[Tl it is enough to show that P^ (limsupt^oc{9j ^t) > O) > 0, or equivalently, 
that Ef, (e"'™™P'T-<5^^*)) < 1. Let e > 0. Use first Fatou's Lemma and then 
Theorem [5] to obtain the estimate 

Ef, (^e-'™'^"P'Too(9,^t>j < liminf P^ (e^^^'^'^) < liminf E^^ (e'^^^*)) (11) 



Jdnos Engldnder / Branching diffusions, superdiffusions and random media 323 

The ergodicity of (i^, P^^) implies that the right hand side of (fTTj) is less than 
one, finishing the proof. □ 

Note the intuition behind the last part of the proof: the spine particle visits 
every part of B for arbitrarily large times because it is an ergodic diffusion; this 
forces the process itself to do the same. 



4.4. Proof of Theorem\B (ii) 

We may assume without loss of generality that A S (0, Ac). By standard theory, 
there exists a large enough B* CC D with a smooth boundary so that 

A* := \ciL + f3,B*) e (A, Ac]. 

Further, we can also choose B* large enough so that supp(//) CC B* . 

We first claim that if flo {limf|oo e~^*Zt {B*) = 00}, then (fio) > 0. 
Indeed, if Z is defined as the (£,/3; _B*)-branching process then 



P^d^o)>P^Aliram^e-^*Zt{B*)>0] > liminf e"^ IZJ > 

y ijoo J y tloo 

> P J \im e~^'' {cb*, Zt) >0] , 

where (L + (3 — X*)ip* = in B* and 0* = on dB*. Note it is implicit in the 
definition of Z that particles are killed on the boundary OB*. Since A* > 0, 
Lemma [1] implies that the last term in the last inequality is positive. 

Now let B be any open set with % ^ B CC D. Let p := inf^g^. a;, B) > 0, 
where {p{t,.,dy) : t > 0} is the transition measure for {Y,F). Let < q < p 
and An := {Zn+i{B) > qZn{B*)}. It follows from the law of large numbers 
and the Markov property that on fio, lim„|oo P{An \ Zn, Zi) = 1. Using the 
extended Borel-Cantelli lemma, it easily follows that \im swp^.^ ^ e"^* Zt{B) = 
00 — a.s. on flo. □ 
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5. The Law of Large Numbers for spatial branching processes and 
superprocesses 

In this section we would like to give a brief review on the problem of the Law of 
Large Numbers for spatial branching processes. Here Z will denote the {L, (3, D)- 
branching diffusion and X will denote the (L, f3, a, -D)-superprocess. 



5.1. The question 

Consider Z and assume that it starts with a (finite) discrete measure /x. Let us 
first assume that Y is conservative (i.e. never reaches the cemetery state A) and 
/3 is a positive constant. Since the total mass process \Z\ is a non-spatial Yule's 
process (pure birth process), it is well known that e~^*\Zt\ has an a.s. positive 
hmit Rst—> oo. That is, \Zt\ ^ ef^Nf,. where N,, > a.s. Let p^{t,x, •) denote 
the kernel corresponding to the diffusion Y (or, equivalently, to L). It is natural 
to ask then whether it is true that for B CC D, 

Zt{B)r^\Zt\-p'^{t,x,B)? 

This would be a kind of Law of Large Numbers, because what it says is that 
the proportion of particles in B at t is given by the probability that a single 
particle is in B at t. In other words, the question is whether 

Zt{B)^e'^*N^-p''{t,x,B), 

and since e^* ■ p^{t, x, •) = p^~^^{t, x, ■) is the kernel corresponding to i + /?, we 
can rewrite this in the form 

Z,{B)^N^-p'^+l^{t,x,B). 

Of course this last formula makes sense even for spatially dependent /3's. 
Moreover, it is known (and easy to derive) that p^+^(i, x, B) is just the expec- 
tation of Zt{B) starting from ji = 5x ('first moment formula'). So, for a general 
/3 one has the following conjecture: 

This conjecture is wrong! The reason is the possibility of local extinction. If 
local extinction occurs, then the left hand side becomes zero in finite time (and 
stays zero). 

In terms of our spectral criterion for local extinction this means that there is 

no Law of Large Numbers, or more generally, there is no scaling for Zt{B) when 
Ac = Ac(I/ + /?) < 0. The correct question therefore is this: is there a dichotomy 
in the sense that either 

(i) Z suffers local extinction (Ac < 0), or 

(ii) the (local) Law of Large Numbers holds for Z(Xc > 0)? 

This is not an easy question! In fact, since the spectral criterion of local 
extinction for superprocesses is the same, one can ask if the same dichotomy 
holds for the superprocess X. Moreover one has to clarify of course, what exactly 
means in the formulas. 
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5.2. Some answers 

There are only partial answers in this direction. For superprocesses, it has been 
proven in [5] that ii D — M.'^ and Ac > and furthermore a certain spectral 
condition holds for L + f3, then for all 7^ g G C^, the ratio ^'\xfg) ^ 
(non-degenerate) limit (depending only on fL) in law. 

One equivalent way of formulating the spectral condition is as follows. In gen- 
eral p'"^^{t, X, B) ^ e"^"^* X possible subexponential term. The spectral condition 
is then equivalent to not having a subexponential term (the scaling is 'purely 
exponential'). For example, when Y is d-dimensional Brownian motion and (3 is 
positive constant, one has p^'^^ {t, x, B) ~ e'^* • t~'^/'^, and therefore this case is 
not included in the setting. On the other hand, if one replaces Brownian motion 
with the Ornstein-Uhlenbeck process (or, in fact, with any positive recurrent 
process), then it is included. 

It is important to point out however that, for nonconstant /3's, the underlying 
motion corresponding to L does not have to he positive recurrent (it can even be 
transient), it is the operator L + P, involving branching too, that has to satisfy 
a certain spectral condition. In fact, an auxiliary diffusion, corresponding to a 
particular /i-transform of the operator L + 3 — Ac is positive recurrent. This 
other diffusion also appears in a 'backboneri construction. 

This result was improved in [B] in the sense that a general Euclidean domain 
£> C M'' was considered instead of Z? = M'' and convergence in law was upgraded 
to convergence in probability; recently it was further extended [4] to certain 
cases when the scaling is not purely exponential, including the case when Y is 
d-dimensional Brownian motion and f3 is positive constant. 

Convergence in probability is far from being satisfying. In fact, by a well 
known theorem T, '5] , the analogous result for branching Brownian motion (with 
constant /3) holds even in the strong sense (i.e. the convergence is a.s.), thus there 
is no real reason to assume that the convergence only holds in probability for 
superprocesses. 

In three very recent papers [21 IS [7] a.s. convergence was finally settled for 
certain classes of branching diffusions [31 and superprocesses P] . Despite the 
somewhat restrictive conditions, [2] seems to be the first time when SLLN is 
demonstrated for superprocesses, and it represents a major step in the direction 
of establishing the 'local extinction vs. SLLN' dichotomy for superprocesses in 
terms of the associated generalized principal eigenvalue. Another nice feature 
of the paper is that besides symmetric diffusions, symmetric Levy processes are 
treated too. An essential difficulty is compounded in the fact that the method 
of [3] for replacing discrete times with continuous times does not carry through. 
The authors of 2 manage to overcome this difficulty with a highly nontrivial 
application of the martingale formulation for superprocesses. 

At this point it seems to be a serious challenge to prove/ disprove the full 
dichotomy (without additional scaling assumptions) for both kinds of processes. 



'Backbone' constructions for superprocesses are similar to the 'spine' constructions for 
branching diffusions. 
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6. Further topics in superprocesses: compact support property and 
polar decomposition 

6.1. The compact support property 

In this subsection we will consider superprocesses with a somewhat more general 
branching mechanism. For simplicity, we take D = M.'^, although most of the 
results are true for general Euclidean domains. 

In this case the evolution equation appearing in the log-Laplace formulation 
is of the form 

ut — Lu + (iu~ auP in R'^ x (0, oo); 

u{x,0)^f{x) inR'^; (12) 

t)>0 in R'' X [0, CO), 

where p G (1,2]. When 1 < p < 2 the variance is no longer finite and only 
the moments up to p exist in the particle process approximation, whereas the 
scaling of the clock becomes nP~^. Finally, the paths of the superprocess are not 
continuous but only cadlag when 1 < p < 2. 

Let /i e Alc(R)- Since D = R'', the measure-valued process corresponding to 
possesses the compact support property if 

Pf, y supp X{s) is bounded = 1, for all t > 0. (13) 
\a<s<t J 

The definition is independent of the choice of the nonzero compactly supported 
starting measure. 

There are four objects, corresponding to four different underlying probabilis- 
tic effects, which can influence the compact support property: 

1. L, the operator corresponding to the underlying motion; 

2. f3, the mass creation parameter of the branching mechanism; 

3. a, the nonlinear component of the branching mechanism, which can be 
thought of as the variance parameter if p = 2; 

4. p, the power of the nonlinearity, which is the scaling power and is con- 
nected to the fractional moments of the offspring distribution in the par- 
ticle process approximation to the measure-valued process. 

We shall see that both L and a play a large role in determining whether or not 
the compact support property holds; /3 and p play only a minor role. 

The connection with partial differential equations is that the compact support 
property turns out to be equivalent to a uniqueness property for solutions to 

m- 

Result 1. The compact support property holds for one, or equivalently all, 
nonzero, compactly supported initial measures fi if and only if there are no non- 
trivial solutions to (|12|) with initial data / = 0. 



The results of this subsection are joint with R. Pinsky. 
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Remark 3. Result [T] suggests a parallel between the compact support prop- 
erty for measure-valued processes and the non-explosion property for diffusion 
processes. Indeed, the non-explosion property for the diffusion process Y cor- 
responding to the operator L is equivalent to the nonexistence of nontrivial, 
bounded positive solutions to the linear Cauchy problem with initial data: 

Ut = Lu in M'' X (0,oo); 

u{x,Q) = in R'^; (14) 
■u > in R'' X (0,oo). 

It is natural for bounded, positive solutions to be the relevant class of solutions 
in the linear case and for positive solutions to be the relevant class of solutions 
in the semi-linear case. Indeed, by Ito's formula, the probabilities for certain 
events related to Y are obtained as bounded, positive solutions to the linear 
equation, and by the log-Laplace equation, the negative of the logarithm of the 
probability of certain events related to X can be obtained as positive solutions 
to the semi-linear equation. o 

The class of operators L satisfying the following assumption will play an 
important role. 

Assumption 1. For some C > 0, 

2. \b{x)\ <C{l + \x\), xeM.'^. 
We have the following basic result. 
Result 2. Let p ^ {^I'A '^'^'^ coefficients of L satisfy Assumption]^ 

1. There is no nontrivial solution to (|14p .' thus, the diffusion process Y does 
not explode. 

2. If 

inf a{x) > 0, 

then there is no nontrivial solution to (|12p with initial data f — 0; thus, 
the compact support property holds for X . 

The conditions in Assumption[T]are classical conditions which arise frequently 
in the theory of diffusion processes. Result [2] shows that if the coefficients of L 
obey this condition and if the branching coefficient a is bounded away from zero, 
then everything is well behaved — neither can the underlying diffusion process 
explode nor can the measure- valued process fail to possess the compact support 
property. 

The following result shows that the compact support property can fail if 
infj-gRd a{x) — 0. It also demonstrates that the effect of a on the compact 
support property cannot be studied in isolation, but in fact depends on the 
underlying diffusion. 



Jdnos Engldnder / Branching diffusions, superdiffusions and random media 



330 



Theorem 4. Let p ^ (li2] and let 



2 dxidxj 

z,j=l ■' 

d 

some Co > 0. 

a{x)>CieM~C2\x\^-^), 

for some Ci, C2 > 0, then the compact support property holds for X . 
2. If 

a{x) < Cexp(-|a;p-'"+^) and (3{x) > -C(l + 

for some C, e > anc? some 5 < e, then the compact support property does 
not hold for X . 

(By Result [T] to prove this theorem, it is necessary and sufficient to show 
that if a is as in part (f) of the theorem, then there is no nontrivial solution 
to (fT2|) with initial data f — 0, while if a is as in part (2) of the theorem, then 
there is such a nontrivial solution.) 

As a complement to Theorem 21 we note the following result. 

Result 3. Letpe (1,2]. 

1. Let d > 2 and let 

L = A{x)A, where A{x) > C(l + l^l)™, for some C > and m> 2. 
Assume that 

sup a{x) < 00 and /3 > 0. 

Then the compact support property does not hold for X . 

2. Let d = I and let 



L = A{x)^j^, where A{x) > C{l + \x\)"^, for some C > and m > 1+p. 



dx 

Assume that 

sup a{x) < 00 and /? > 



Then the compact support property does not hold for X . 
3. Let d — 1 and let 



da 

Assume that 



L^A{x)-—r, where A{x) < C{l + \x\)"^, for some C > 0, andm<l+p. 
dx 



inf a{x) > and j3 <0. 
Then the compact support property holds for X . 
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Remark 4. It follows from Result [3] that if d 1 and L = (1 + |a;|)™£7, with 
m g (2,3], and say a = 1 and P — 0, then the compact support property will 
depend on the particular choice of p e (1, 2]. o 

Remark 5. If d = 2 and L = (1 + |x|)™A, with m > 2, and say a = 1 and 
/3 = 0, then by Result [U the compact support property does not hold, yet the 
underlying diffusion does not explode since it is a time-change of a recurrent 
process; namely, of two-dimensional Brownian motion. o 

The proof is based on the fact that, by Result [1] it is necessary and sufficient 
to show that under the conditions of parts (1) and (2), there exists a nontrivial 
solution to while under the conditions of part (3) there does not. 

A useful comparison result is as follows. 

Proposition 3. Assume that 

/3i </32 

and 

< a2 < Qfi. 

// the compact support property holds for (3 — P2 o-nd a — ai, then it also holds 
for P — Pi and a = ai. 

Theorem m and Result [3] demonstrate the effect of the underlying diffusion Y 
on the compact support property in the case that L is comparable to (H-I^D^A. 
We now consider more generally the effect of the underlying diffusion process 
on the compact support property. 

Theorem 5. Let p G (1,2] and assume that the underlying diffusion process Y 
explodes. Assume that 

sup a{x) < 00 and inf P{x) > —00. (15) 
Then the compact support property does not hold. 

Remark 6. Theorem [5] shows that if the branching mechanism satisfies 
then the compact support property never holds if the underlying diffusion is 
explosive. The converse is not true — an example was given in Remark[2]following 
Result [21 and another one appears in the remark following Corollary [31 o 

The next result shows that the restriction sup^gjjd a{x) < cxd in Theorem [5] 
is essential. 

Proposition 4. Let pE [1,2]. Let m E {—00,00), 

i = (l + |a;|)™A in R'^, 

P = and 

a{x) > c{l + \x\)"'~\ 

for some c > 0. Then the compact support property holds for the measure-valued 
process X. However, if m > 2 and d>3, the diffusion process Y explodes. 
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We now discuss how the concept of a ineasure-vahied process hitting a point 
can be formulated and understood in terms of the compact support property. 

Similarly to the range define also TZt = { U^gjo^tj supp(X(s))). A path of the 
measure-valued process is said to hit a point xq G K'' if xq G TZ. li X exhibits 
local extinction, then G 7?. if and only if xq G TZt for sufficiently large t. Thus, 
we have: 

If X exhibits local extinction, then 

Pfj,{X hits a;o) > if and only if there exists a, t > such that (16) 
Pfj.i^o<s<t supp(X(s)) is not compactly embedded in R'* — {xq}) > 0. 

Now although we have assumed in this section that the underlying state space is 
Mf^, everything goes through just as well on an arbitrary domain D C M.'^. Recall 
that the compact support property is defined with respect to the domain D, and 
the underlying diffusion will explode if it hits A in finite time. In particular, 
Resuh [U still holds with R'^ replaced by D. 

In light of the above observations, consider a measure-valued process X cor- 
responding to the log-Laplace equation on M.'^ with d > 2. The underlying 
diffusion process Y on M'* corresponds to the operator L on R'^. Let Y denote 
the diffusion process on the domain D = M.'^ — {xq} with absorption at xq and 
corresponding to the same operator L. (Note that if xq is polar for Y, then 
Y and Y coincide when started from a; 7^ xg. In fact, xq is always polar un- 
der the assumptions we have placed on the coefficients of L.) Let X denote 
the measure- valued process corresponding to the log-Laplace equation (|12p . but 
with M'' replaced by Z? = R'' — {xq}. It follows from (fT6| that if X suffers local 
extinction, then the measure-valued process X hits the point xq with positive 
probability if and only if X on — {xq} does not possess the compact support 
property. Furthermore, the above discussion shows that even if X does not suf- 
fer local extinction with probability one, a sufficient condition for X to hit the 
point Xq with positive probability is that on R"* — {xo} does not possess the 
compact support property. 

A similar analysis can be made when d = 1. (The process X above must 
be replaced by two processes, X+ and X~ , defined respectively on (xo,oo) and 
(— C!0,a;o).) In the sequel, we will assume that d>2. 

Consider now the following semi- linear equation in R'' — {0}: 

ut = ^Au - uP in (M'' - {0}) X (0, 00); 

u{x,0) = in R'^ - {0}; 
u > in (M'' - {0}) X [0,00). 

By Result [TJ the measure-valued process corresponding to the semi-linear 
equation ut — ^Au — in R'' — {0} will possess the compact support property 
if and only if ([TT]) has a nontrivial solution. The following result was recently 
proved by P. Baras, and M. Pierre ^1,, and later by R. Pinsky [B], using a different 
method. 



(17) 
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Result 4. Let p > 1 and d >2. 

1. If d < -^zj: then there exists a nontrivial solution to (|17p . 
2- Ifd>^, then there is no nontrivial solution to (jl7p . 

Here p > 1 is unrestricted, even though of course there is no probabiHstic 
import when p > 2. 

In terms of the superprocess, the foUowing was shown in [4] (see also the 
remark at the end of this subsection). 

Theorem 6. Let p > 1 and d > 2. Let X denote the measure-valued process on 
the punctured space R"^ — {0} corresponding to the semi-linear equation Ut = 
iAu -\- fiu — in (R'^ — {0}) x (0, oo). Consider the Cauchy problem 

Ut = ^Au + (3u-uP in [W^ - {0}) x (0, oo); 

u(a;,0) = m M"* - {0}; (^^) 
u>{) in{«^ - {0}) X [0, oo . 



Assume that 



Let 



d < 



p-1 



1- If 



d{p-l)-2p 



> ^7— fT"-, for some k G (0, — /3o], 



then there exists a nontrivial solution to (|18p ; hence, the compact support 
property does not hold for X . 
2. If 

limsup |a;p/9(x) < /3o, 



then there is no nontrivial solution to (jlSp ; hence, the compact support 
property holds for X . 

Remark 7. The restriction k < —(3o is made to ensure that /? is bounded from 
above. o 

We have the following corollary of Result [3] and Theorem [SI 

Corollary 3. Let X denote the measure-valued process on all of R'^ , d > 2, 
corresponding to the semi-linear equation Ut — ^Am + fiu — in x (0, oo). 

1. If (3 is bounded from below and d < then X hits any point with 

positive probability; 
2- If P <0 and d > , then X hits any point xq with probability 0; 
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3. If P < 0, d < and (3 has a singularity at the origin such that 

limsup|xp/3(x)<fc^l^, 

then X hits with probability 0. 

When P = 0, part (1) follows immediately from Result [H Result [T] and the 
discussion preceding Result 21 The general case requires a little more works as 
one needs to compare different /3's. (Recall that we are always assuming here 
that /? is bounded.) 

When /5 < 0, it is well-known that the super-Brownian motion in the state- 
ment of the corollary becomes extinct. Thus, part (2) follows from Result 21 
Result [H and the discussion preceding Result 21 while part (3) follows from 
Theorem |6l Result [1] and the discussion preceding Result 21 

Remark 8. When /3 = 0, the results in parts (1) and (2) of Corollary O state 
that critical, super-Brownian motion hits a point with positive probability if 
d < -^zi, and with zero probability if d > (Note that this also yields an 

example where X does not possess the compact support property even though 
the underlying diffusion process does not exploded) This result was first proved 
by D. Dawson, I. Iscoe, and E. Perkins i2j, and by E. B. Dynkin and S. Kuznetsov 
[3]. Their methods were different than the one presented here. o 



6.2. Polar decomposition 

As we have already seen in the beginning of the previous section, the total mass 
process of a branching diffusion is sometimes easier to handle than the whole 
process. This is true for superprocesses too. In fact, one can even think of the 
superprocess until its extinction as a process of two components: the total mass 
process, as 'radial component' and the process normalized by the total mass 
(which is a probability measure valued process) as the 'angular component'. 
We are going to demonstrate this for the critical super-Brownian motion (the 
(iA, 0, 1, i?'')-superprocess). 

To this end, let us have a review first on the Fleming- Viot superprocess which 
is of independent interest too. Imagine that we have N individuals on the lat- 
tice Z"* and each of them perform independent symmetric random walk. In the 
genetic setting is the space of genetic types and the random walks are mu- 
tations. Each individual has her own independent exponential clock with rate 
7 > {sampling rate) and when the clock rings, she jumps to a position ran- 
domly chosen from the current empirical distribution of the population [sam- 
pling mechanism). Note that we can consider this mechanism as branching: one 
individual dies and another one is born somewhere. Also, the process can be 
viewed as a (discrete) measure valued process. 



It is easy to find the connection and we leave this to the reader. 



Jdnos Engldnder / Branching diffusions, superdiffusions and random media 



335 



Let us rescale now the model as follows: assign mass to each individual 
(and so the measure valued process becomes a probability measure valued pro- 
cess) . Furthermore speed up time by factor N and reduce lattice mesh by factor 
(So the mutation part will converge to Brownian motion). 

Formally, let VN{t, A) :— J2j/^eAPi'^^'^)' where A C Z'^ and p{s,j) is the 
proportion of the population at the point j at time s > 0. Let N —^ oo now. If 
Viv(0, •) converges to the probability measure v on weakly, then the process 
{VN{t,-);t > 0} converges to a process {V{t,-);t > 0} weakly. This limiting 
(probability measure valued) process is the Fleming- Viot superprocess. 

Let us write X in a 'polar coordinates': Xt = rt ■ dt, where rt := \Xt\ and 
9t j^J-j- for rt > 0, (otherwise rt = 0,Xt — 0). Then it is easy to show 
that r is a diffusion process ('Feller's diffusion) on D = [0,oo) corresponding to 
the operator L = (This diffusion will hit zero with probability one and 

afterwards it stays at zero.) A result due to P. March and A. Etheridge then 
states that is & Fleming- Viot superprocess, except that the sampling rate is 
replaced by l/rt- This 'skew-product' representation is the polar decomposition 
for critical super-Brownian motion. 
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7. Classical problems for random media: Brownian motion among 
Poissonian obstacles and related topics 

In this section we give a brief (and surely not comprehensive) review on some 
of the classical problems in random media. 

7.1. Brownian survival among Poissonian traps 

As usual, B{x, r) denotes the open ball of radius r around x. Let W denote the 
d-dimensional Brownian motion (Wiener-process) starting at the origin (the 
corresponding probability will be denoted by P) and let us define the Wiener- 
sausage up to time t with radius a: 

V,^ = Vo\i y BiWua)] . 

\a<s<t J 

It turns out that this object is intimately related to the problem of Brownian 
survival among Poissonian traps. Let u> denote the Poisson point process (PPP) 
with intensity v{dx) on . (The corresponding probability will be denoted by 
P.) That is, the number of points in ^ C M'^ bounded is Poisson with parameter 
y{A) and the number of points in A c M*^ and the number of points in B c 
are independent random variables if A and B are disjoint. So for example 
V{A is empty) = e~'''^^\ Unless indicated otherwise we will set v{dx) = v ■ 
dx, v > 0. 

The (random) set 

K:= y B{xi,a) 

XiGsupp(aj) 

is a trap configuration, or obstacle. 

If T denotes the first hitting time of K by our Brownian motion, then we 
would like to know what the asymptotic behavior of the probability of {T > t} is. 
Of course this question can be asked in two different ways: what the asymptotics 
of (E (g) P)(T > t) is, or what the P-a.s. asymptotics of P{T > t) is. In the first 
case one talks about the annealed asymptotics and in the second one, about the 
quenched one. 

Now, the annealed asymptotics is in fact just the same as the asymptotics 
for the Laplace transform of the Wiener sausage. Indeed, 

(E O P)(r > t) = (E O P) I y B{Wt,a) is empty = Ee-''^^\ 

\0<s<t I 

A little more than 30 years ago Donsker and Varadhan proved that for all 
a, 1/ > 0, 

lim t-'^/('^+2)logi;exp(-z/14'^) = -~c{d,v), 
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where c(d, v) is an explicitly given positive constant (we are using c because 
there is another constant c{d,i') in the quenched case). Equivalently, 

Ee-"^^" = exp[-£(d, iy)t'^^^''+^\l + o(l))]. 

Note that the scaling is milder than exponential, it is often referred to as 
^stretched exponentiaP. For example in two dimension the Laplace transform 
behaves like 6-==°""'^. 

Of course, the same is then true for the annealed survival probability: 

(E (g,P){T>t)^ exp[-5(d, z.)i'^/('^+2)(l + o(l))]. 

The easy part is 

liminft-''/(''+2) \ogEexp{-vVt'') > -c{d,iy), 

t — »oo 

because all one has to do is to 'advise' a survival strategy which is realized (at 
least) with probability exp[— c((i, z^)i'^/'-'*+^^ (1 + o(l))], whereas in order to verify 

limsupt"'^/('*+2) \ogEcyiY>{-vVt^) < -c{d,iy), 

t — *oo 

one has to show that the advised strategy is actually optimal. 

In the context of Brownian survival among Poissonian traps these (and also 
the quenched asymptotics) were settled by A. S. Sznitman. 

Let us see now how to 'advise' a survival strategy. An obvious observation is 
that if -6(0, Rt + a) is empty (such an empty ball is called a clearing), and if W 
is confined to i?(0, Rt) up to t, then T > t. That is, 

(E (g)P)(T>t)> P(B(0, Rt + a) is empty) • P(rB(o,i^,) > t), 

where Ts(o,i?t) denotes the exit time of W from the ball. It is not hard to 
show that P{TB(o,Rt) > t) ^ c{d)e~^'^^^^'* , where Xd is the principal Dirichlet 
eigenvalue of —A on the unit ball. Hence 

(E ®P){T>t)> c(d)e-'^"''(^'+")'-^''/-«'-*, 

where LOd is the volume of the d-dimensional unit ball. So 

log(E ® P){T >t)> logc(d) - vLOdiRt + a)^ - Xd/Rt ■ t. (19) 

Of course, one wants to get the optimal (strongest) lower estimate. In order 
to achieve this, we match the orders of the last two terms (the constant will 
be lower order and thus negligible), which happens when Rt = C'(t^/('^+^)). 
If Rt = fci^/('*+^', then we have to maximize f{k) :— vujdk'^ + which is 

/ „, \ 1/(^+2) 

found easily (by differentiation): = ( ^/jj^ 1 . The rest is straightforward. 

Setting Rt := ( J^j • 1^/'^'^+^'^ in dividing by t'^/id+'^) and finally 

taking liminf gives the desired lower estimate. 
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As we have mentioned, the proof of the upper estimate is much harder and 
it is definitely beyond the scope of these notes. 

As far as the quenched problem is concerned, the asymptotics is different. If 
a is sufficiently small such that the origin belongs to an infinite trap free cluster 
with positive P-probability, then on that set P-a.s, 

P{T >t)^ exp (^-c(d, (^^^(1 + «(1))) , 

where again c{d, v) is an explicitly given positive constant. 

Finally let us note that the annealed and quenched asymptotics hold also for 
soft obstacles (so far we have been talking about hard ones) , where the meaning 
of being 'soft' is as follows. Let 6 > (the size of b plays no role) and for a given 
realization of u let 

P{T >t)^ exp (^-bj^ txiW,) ds^ . 

(This is not the usual formulation because we do not add up the rates on over- 
lapping balls, but it does not matter from the point of view of asymptotics 
and this is what we will need.) That is, as long as W is inside K it 'feels' an 
exponential killing rate, whereas in the 'Swiss cheese' K'^ it is 'safe'. 

7.2. Random walk in random environment 

Let us first assume d = 1. There are two basic models for random walk in 
random environment (RWRE) : 

• Site randomness (this is what usually meant by RWRE). For a given 
realization of the environment w, {^n} is a Markov-chain on Z. The walker 
located at a; e Z moves to a; -t- 1 with probability p{x, lu) and to a; — 1 with 
probability q{x,Lu) — 1 -~p{x,uj), where {p{x,uj), a; G Z} are i.i.d. random 
variables with values in [0, 1]. 

• Bond randomness (sometimes called 'random conductivity' model). Here 
{cx,x+i{^), a: e Z} are i.i.d. with values in (0, oo) and for a given realiza- 
tion {^n} is a Markov-chain on Z with probability 

/ \ Cx^x+l{l^) 

P(x,uj) := -— -- 

Cx^l,x{^) + Cx,x+l[t^) 

of jumping to a; -I- 1 and with probability q{x,uj) — 1 —p{x,uj) of jumping 
to a; — 1, in which case the {p{x, lu), a; € Z} are of course no longer i.i.d.. 

Site randomness goes back to Chernov |2j in the early sixties, while bond ran- 
domness goes back to Fatt [3] in the fifties. 

Of course it is no problem to define these models in higher dimensions 
too. However, the higher dimensional case is significantly harder than the one- 
dimensional one, and much less is known. In the sequel we keep assuming d = 1. 
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One of the natural questions concern the speed of the random walk. Now, in 
the bond randomness model one can show that for almost all realizations of w, 
^ tends to zero, which is interpreted as 'vanishing limiting velocity'. However, 
for the site randomness case, the situation is far more subtle. Here we suffice 
with citing the classical result of Solomon [4j . 

Theorem 7. Let p{x,uj) :— and p :— p(0,w). // P corresponds to the 

environment and to the walk given lu, then there are three cases. 

1. E[log/9] < and lim„X„ = +oo, 

2. E[log/9] > and lim„X„ = — oo, 

3. E[log/9] — and liminf„X„ = — oo hut limsup„X„ = +oo, 

where all the limits are meant P (g) P^-a.s. 

Furthermore, — ^ v V ® P^-a.s. as n — > oo, where, again, there are three 
cases. 

1. Ep < 1 and V ^ e (0, 1), 

2. Ep > 1, E[p-i] >l andv = Q, 

3. Ep-i <1 andv^ ^ (-1,0). 

Remark 9. To understand the statements better, note the following. 

(i) In the third case Ep~^ < 1 automatically implies Ep > 1, by Jensen. Also 
by Jensen, Ep < 1 implies E[logp] < and Ep~^ < 1 implies E[logp] > 0. 

(ii) When u = 0, it is known that the propagation can be as slow as (logn)^ 
per n steps. 

(iii) It is easy to show (again by Jensen), that Ep < 1 implies w < Ep — Eg, 
where p — p{0, lu) and q — 1—p. Similarly Ep^^ < 1 implies v > Ep — Eg. These 
inequalities are strict (unless p is deterministic), which shows the 'slowdown' of 
the walk relative to the naive guess that the speed is simply v* :— Ep — Eg. o 
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8. Spatial branching processes among Poissonian obstacles: hard 
obstacles 

We now start putting Poissonian obstacles and spatial branching togetherP^ 

Let T denote the hitting time of K (defined in the previous section) by 
Z, where Z denotes branching Brownian motion (BBM) with strictly dyadic 
branching and constant rate f3 > 0, starting with a single particle at the origin. 
That is T is the time the first particle hits any trap: 

T := inf{s > 0, Ws e K}. 

Our question is again the asymptotic behavior of the probability of {T > t}. 

8.1. Constant trap intensity 

Let P be as in the previous section and let P denote the probability correspond- 
ing to the BBM. 

Theorem 8. When d > 2, the annealed asymptotics is as follows: 
lim - log(P (E)P)(T>t) = -(3, 

t — too t 

and this is also true for the supercritical super- Brownian motion starting with 
unit mass at the origin. More precisely, if X is the {^A, l3,a;R'^)-superprocess 
with a, (3 positive constants and d > 2, then 

lim i log(P ® Ps,){T* >t\S)^-P, 

t— too t 

where P corresponds to X, T* :— inf{t > | Xt{K) > 0} and S is the event 
that X survives forever. 

Although this result is far from being trivial, in a sense it is not too exciting 
because it indicates the overwhelming influence of the branching over the spatial 
motion. For example, in case of Z , the probability of not branching at all up to t 
is e~^* and by the Donsker-Varadhan result we know that a single particle avoids 
K up to t with a subexponentially small probability. Therefore, we immediately 
obtain the lower bound 

lim inf - log(P (g) P){T > t) > -/?, 

t— too t 

for all d > 1. That this actually coincides with the upper bound for d > 2, 
shows that in a sense suppressing the branching is the 'best the system can do' 
in order to avoid the traps, and the cost of doing so is the dominating factor 
on a logarithmic scale. Because of the conditioning for X however, the lower 
estimate for the superprocess is not easy. 



^^Results of the first subsection are due to the author, results in the second one are joint 
with F. den Hollander. 
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8.2. Decaying trap intensity 

In order to treat the more delicate one-dimensional case of the previous sub- 
section, we are going to discuss a setting when the trap intensity decays at the 
order |a;|''~-'^ as |x| oo, which includes arc previous one- dimensional case. So 
let us consider now a trap intensity satisfying that 

where ^ is a fine tuning constant. Considering the same process Z under P as 
before and the same killing rule (we kill Z when the first particle hit K) we 
have the following annealed asymptotics for the tail of the survival probability. 

Theorem 9. 

lim logi(P® P)(r >t) = -I{e,(3,d), 

t — *oo t 

where T is as before and (3,d) > can he computed in terms of a variational 
problem. 

Instead of writing down the variational problem, let us see the optimal sur- 
vival strategy the variational problem 'encodes'. 

First, it turns out that there is a crossover at a critical value of £, which 
we will denote by (.cr = ^cr{P,d). (This crossover is the reason we chose this 
particular order for If the density is higher or lower order than |a;|^~'^, then 
we fall into one of the regimes already present when ^ ~ £|a;|^^'', and i > £cr 
OT £ < icr, respectively.) If := s^^-\//3/2, where Sd is the surface of the d- 
dimensional unit ball, then for d= 1, l^r = (*cri while for d>2, Icr < (*cr- the 
theorem below we try to explain the optimal survival strategy in words rather 
than by using the precise formulation. Roughly speaking, the optimal survival 
strategy means that, conditioned on survival, the system, consisting of the BBM 
as well as the traps, behaves in a certain way up to time t with P(g)P-probability 
tending to one as t ^ oo. We will distinguish between a 'low intensity' {£ < £„) 
and a 'high intensity' {£ > £cr) regime. Note, that the case £ = £cr is left open. 

Concerning the optimal survival strategy, we have the following result. (The 
first part is actually just an interpretation of the precise statement which will 
be given in the theorem after this result.) 

Result 5. Conditional on {T > t}, the following holds. 

• (Low intensity regime) For £ < £cr, the ball B{0, \/2pt) is emptied, the 

system stays inside this ball and, branches at rate (3. 

• (High intensity regime) For £ > £cr, 

— d= 1: a B{0, o{t)) ball is emptied, while Z suppresses branching until 
t and stays inside this ball, 

— d > 2: There are two numbers (these are actually the unique min- 
imizers of the variational problem for I{£, P,d)), < rj* < 1 and 
c* > 0, such that Z suppresses branching until rft, empties a ball of 
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radius — ri*)t around a point at distance c*t from the origin; 

during remaining time {1 — rf)t the BBM branches at its normal (3 

rate. 

Furthermore, for d> 2, £ i-^ I{£, /?, d) is continuous and strictly increasing, with 
lime^^ I {i,P,d) = (3 (see Fig. 1) and I ^ ri*{£,P,d) and £ c*{£,j3,d) are 
both discontinuous at £cr and continuous on {icr, oo), with 

lim =1, lim c*(£,P,d) = 0. (20) 

e^oo c*{l,(3,d) ' e->oo ^ ^ ' 

Finally, c* > - rf) for all l> £„■ 




Fig. 1 £^ I{t,(3,d) for: (i) d = 1; (ii) d>2. 



Remark 10. By the last statement, the clearing never contains the origin for 
£ > £cr- However, for £ < £cr, it does. So, the discontinuity mentioned in the 
theorem means that, surprisingly, when £ crosses its critical value, both c* (the 
distance of the clearing normalized by t) and \/2/3(l — r?*) (the radius of the 
clearing normalized by t) have a jump. So in a sense, the clearing 'jumps away' 
from the origin and its size 'jumps down'. We have no intuitive explanation for 
this phenomenon. o 

The precise statement regarding the optimal survival strategy is as follows. 

Theorem 10 (optimal survival strategy). Fix l3,a. For r,b > and t > 0, 

define 

C{t;r,b) = {3x0 eW^: \xo\ = b, B^tixot) H K = 9}. (21) 
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(i) For d= 1, i < £cr or d>2, any I, and < e < 1 - 

^lim (E X Ps,) (C (t; ^ (1 - r;* - e) , c*) | T > = 1, 
lim (E X PaJ f \Z(t)\ > [e^i^-ri'-^)*^ \ T>t]=l. (22) 

(ii) For d> 1, i < £cr and e > 0, 

^lim (E X P,j(B(^^^,^^(0) ni^ ^ I T > = 1, 

^lim (E X P,j(i?(i) C B(,^,)^,(0) \T>t)=l, 

^lim (E X Ps,)[R{t) ^ B(^_^^^^(0) \T>t)=l. (23) 

frnj For d > 1, £ > 4r and < e < t]* , 

lim (E X PsJ\Z{{v* - e)t)\ < [t'^+^J \T>t)=l. (24) 

(iv) For d = 1, £ > £cr and e > 0, 

lim (E X Ps,)(B,t{0) nK ^ID\T>t) =1, 

lim (E X P5o)(R{t) C B,t{0) \T>t)=l. (25) 

When d — 1 wc get an answer for the one dimensional constant intensity 
case. In fact, the analysis shows that the optimal strategy is more extreme: in 
the low intensity regime I{£, /3, d) grows linearly with £ and reaches the value /?, 
while in the high intensity regime it is always /3 (see Fig. 1). This corresponds 
to the picture that in the low intensity regime there is neither suppressing of 
branching nor dislocation (ry* = c* =0), but an interval around the origin must 
be emptied, whereas in the high intensity regime, the branching is completely 
suppressed (77* = 1, c* =0). 

The following problems are open: 

- For d = 1 and £ > £cr, what is the radius of the o(t)-ball that is emptied 
and how many particles are there inside this ball at time t? 

- For d > 2 and £ > £cr, what is the shape of the "small tube" in which 
the system moves its particles away from the origin while suppressing the 
branching? How many particles are alive at time r/*t? 

- What can be said about the optimal survival strategy at £ = Icr^ 

- Instead of letting the trap density decay to zero at infinity, another way 
to make survival easier is by giving the Brownian motion an inward drift 
while keeping the trap density field constant. Suppose that dv/dx = £ and 
that the inward drift radially increases like K;|a;|''~^, |a;| ^ 00, k > 0. Is 
there again a crossover in £ at some critical value £cr = £cr{i^, d)l What 
is the optimal survival strategy? 
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We close this section with discussing briefly another killing rule, which we 
call 'individual killing rule'. We now only kill Z when all its particles have been 
absorbed. So the death time becomes 

f :=inf{i>0| |Zf|=0}, 

where is the BBM with killing at dK. 

What can we say about the tail behavior of the survival probability then? 
Well, actually there is no tail behavior because there is no tail at all! And this 
is true for any locally finite trap intensity! 

More precisely, the following holds. 

Theorem 11. Fix d,(3,a. For any locally finite intensity measure v, 

lim (E X P5j(f > t) > 0. (26) 

t — >oo 

Theorem [Tl] follows from the assertion that the system may survive by empty- 
ing a ball with a finite radius R> Rq, where Rq is chosen such that the branching 
rate /? balances against the killing rate X{B{0, Rq)), the principal Dirichlet eigen- 
value of —A/2 on the ball B{0, Rq): X{B{0, Rq)) — (3. Indeed, it can be shown, 
that for any R > Rq there is a strictly positive probability (denoted by pn) that 
at all times at least one particle has not yet left -6(0, Rq). Consequently, the sur- 
vival probability is bounded from below by sup^^j^^ {pji exp[— z^(i3(0, i?o))]}• 
The following two problems for individual killing are open: 

- What is the limit in ((26)) . say, when diz/dx is spherically symmetric? 

- If the Brownian motion is given an outward drift, then for what values of 
the drift does the survival probability decay to zero? 
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9. Spatial branching processes among Poissonian obstacles: 'mild' 
obstacles 

9.1. Basic notions and questions 

Let us consider now a mechanism which is 'milder' than kiUing but still reduces 
the population size. The purpose is to study a spatial branching model with 
the property that the branching rate is decreased in a certain random region. 
More specifically we will use a natural model for the random environment: mild 
Poissonian obstacles. 

We arc going to work with slightly different notation relative to the previous 
sections. As before, let w be a Poisson point process (PPP) on M'' with intensity 
v > and let now P denote the corresponding law. Furthermore, let a > and 
< /3i < /32 be fixed. We define the branching Brownian motion (BBM) with 
a mild Poissonian obstacle, or the '(z^, /32, a)-BBM' as follows. Let B{z,r) 
denote the open ball centered at z G R'^ with radius r and let K denote the 
random set given by the a-neighborhood of w. 

K = K^:= U B{xi,a). 

XiGsupp(w) 

Then K is a mild obstacle configuration attached to uj. This means that given 
u, we define P'^ as the law of the strictly dyadic (i.e. precisely two offspring) 
BBM on K*^, d> 1 with spatially dependent branching rate 

l3{x,oj) := PilK^ix) + P2lKs,{x). 

An equivalent (informal) definition is that as long as a particle is in K'', it obeys 
the branching rule with rate P2, while in K its reproduction is suppressed and 
it branches with the smaller rate (We assume that the process starts with 
a single particle at the origin.) The process under P'^ is called a BBM with 
mild Poissonian obstacles and denoted by Z. The total mass process will be 
denoted by \Z\. Further, W will denote d-dimensional Brownian motion with 
probabilities {P^,, xG R''}. 

Remark 11 (self-duality). The discrete setting has the advantage that when 
is replaced by Z'', the difference between the sets K and K"^ is no longer 
relevant. Indeed the equivalent of a Poisson trap configuration is an i.i.d. trap 
configuration on the lattice, and then its complement is also i.i.d. (with a dif- 
ferent parameter). So, in the discrete case 'Poissonian mild obstacles' give the 
same type of model as 'Poissonian mild catalysts' would. This nice duality is 
lost in the continuous setting as the 'Swiss cheese' K'^ is not the same type of 
geometric object as K. o 

Consider now the following natural questions (both in the annealed and the 
quenched sense): 

1. What can we say about the growth of the total population size? 
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2. What are the large deviations? (E.g., what is the probability of producing 
an atypically small population.) 

3. What can we say about the local population growth? 

As far as (1) is concerned, recall that the total population of an ordinary 
(free) BBM grows a.s. and in expectation as e^^*. (Indeed, for ordinary BBM, 
the spatial component plays no role, and hence the total mass is just a /32-rate 
pure birth process X. As is well known, the limit N := limj^oo e~^'^*Xt exists 
a.s and in mean, and P(0 < N < oo) = 1.) In our model of BBM with the 
reproduction blocking mechanism, how much will the suppressed branching in 
K slow the global reproduction down? Will it actually change the exponent /32? 
(We will see that although the global reproduction does slow down, the slowdown 
is captured by a sub-exponential factor, being different for the quenched and 
the annealed case.) 

Consider now (2). Here is an argument to give a further motivation. Let us 
assume for simplicity that (3i = and ask the simplest question: what is the 
probability that there is no branching at all up to time t > 01 In order to 
avoid branching the first particle has to 'resist' the branching rate /32 inside 
K'^. Therefore this question is quite similar to the survival asymptotics for a 
single Brownian motion among 'soft obstacles' — but of course in order to 
prevent branching the particle seeks for large islands covered by K rather then 
the usual 'clearings'. In other words, the particle now prefers to avoid the if" 
instead of K. Hence, (2) above is a possible generalization of this (modified) 
soft obstacle problem for a single particle, and the presence of branching creates 
new type of challenges. 

Finally, for further mathematical models in population biology see e.g. [4]. 



9.2. Some results on the population growth 



Let us consider the total population size first. Since the behavior of the total 
population size does not seem to be easy to handle, let us first try to handle its 
expectation and then trusting in the Law of Large Numbers (i.e. that the total 
population size behaves as its expectation) , we will at least know, what we want 
to prove. 

Concerning the expected global growth rate we have the following result, 
which, although not at all trivial, can very easily be derived from the Donsker- 
Varadhan result. 



Lemma 2 (Expected global growth rate). On a set of full V -measure, 

t 



exp 



P2t — c{d, v)- 



Ki + o(i)) 



, as t 



ilogtf/d 
(quenched asymptotics), and 

{^®E^)\Zt\^eMP2t-~c{d,v)t'^'^'^^''\l + o{l))], ast 



(27) 



(28) 
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(annealed asymptotics) , where 




) 



-2/d 



c{d, v) := \d 

\vujd 



( 



d + 2 
2 



) 



( 



2Ad 
d 



) 



d/(d+2) 



c((i, v) 



and ujd is the volume of the d- dimensional unit ball, while Xd is the principal 
Dirichlet eigenvalue o/— on it. 

Notice that 

1. (3i does not appear in the formulas, 

2. the higher the dimension is, the smaller the expected population size is. 

Remark 12. Let us pretend for a moment that we are talking about an ordinary 
BBM. Then at time t one has e^^* particles with probability tending to one as 
t — > cx) (the population size divided by e''^* has a limit, to be precise). For t 
fixed take a ball B = B(0, R) (here R = R{t)) and let K be so that B ^ K" 
(such a ball left empty by K is called a clearing) . Consider the expected number 
of particles that are confined to B up to time t. These particles do not feel the 
blocking effect of K, while the other particles may have not been born due to 
it. 

Optimize R{t) with respect to the cost of having such a clearing and the 
probability of confining a single Brownian motion to it. This is precisely the same 
optimization as for the classical Wiener-sausage. Hence one gets the expectation 
in the theorem as a lower estimate. 

One suspects that the main contribution in the expectation in (j28|) is coming 
from the expectation on the event of having a clearing with optimal radius R{t). 
In other words, denoting by pt the probability that a single Brownian particle 
stays in the i?(i)-ball up to time t, one argues heuristically that ptc^^*" particles 
will stay inside the clearing up to time t 'for free' (i.e. with probability tending 
to one as t — » oo). 

The intuitive reasoning is as follows. If we had independent particles instead 
of BBM, then, by a 'Law of Large Numbers type argument' (using Chebysev 
inequality and the fact that limt_>oo Ptc'^^ — cxd), roughly pte'^^* particles out of 
the total e'^^* would stay in the i?(i)-ball up to time t with probability tending 
to 1 as t t OO- One suspects then that the lower estimate remains valid for 
the branching system too, because the particles are 'not too much correlated' . 
This kind of argument (in the quenched case though) can be made precise by 
estimating certain covariances. 

To understand the difference between the annealed and the quenched case 
note that in the latter case large clearings (far away) are automatically (that is, 
P-a.s.) present. Hence, similarly to the single Brownian particle problem, the 
difference between the two asymptotics is due to the fact that even though there 
is an appropriate clearing far away P-a.s., there is one around the origin with 
a small (but not too small) probability. Still, the two cases will have a similar 
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element when one drops the expectations and investigate the process itself, and 
show that inside such a clearing a large population is going to flourish. This 
also leads to the intuitive explanation for the decrease of the population size as 
the dimension increases. The radial drift in dimension d is (d— l)/2. The more 
transient the motion is (that is, the larger d is), the harder it is for the particles 
to stay in the appropriate clearings. o 

Let us prove now the expectation result. First we need a lemma. 

Lemma 3 (Expectation given by Brownian functional). Fix lo. Then 



E'^\Zt\ =Eexp 



ds. (29) 



Proof. It is well known ('first moment formula' of spatial branching pro- 
cesses), that -E^l^tl = {Ttl){x), where u{x,t) := (Ttl)(x) is the minimal solu- 
tion of the parabolic problem: 

^ = Qa + Z^^wouR'^x (0,(X3), 

u{;0) = 1, (30) 
u > 0. 

(Here {Tt}t>o denotes of course the semigroup corresponding to the generator 
+ (3 on R''.) This is equivalent (by the Feynman-Kac formula) to (P^ . □ 
Proof of Lemma [2j: Since (3 := PiIk + h^K" = h — {P2 — Pi)tK, we can 
rewrite the equation (|29|) as 



E'^\Zt\= e^'^'Eexp 



{[32-pi)tK{Ws)ds 



The expectation on the righthand side is precisely the survival probability among 
soft obstacles of 'height' f32 — Pi- (As we have already mentioned, that one 
does not sum the shape functions on the overlapping balls, does not make any 
difference with regard to the asymptotics) . The statements thus follow from 
the Donsker-Varadhan asymptotics for soft obstacles . □ 

We now investigate the behavior of the (quenched) global growth rate. 

As already mentioned in these notes, it is a notoriously hard problem to prove 
the Law of Large Numbers in full generality for spatial branching systems, and 
the not purely exponential case is particularly challenging. 

To elucidate this point, let us consider the (L, (3, D)-branching diffusion with 
some Euclidean domain D C M.'^. Let ^ / be a nonnegative compactly sup- 
ported bounded measurable function on D. If Ac, the generalized principal eigen- 
value of i -I- /3 on £) is positive and S = {St}t>o denotes the semigroup corre- 
sponding to L+P on D, then {Stf){-) grows (pointwise) as e"*^"* on an exponential 
scale. However, in general, the scaling is not precisely exponential due to the 
presence of a subexponential term. 
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In the present case it turns out that the resulting operator + does not 
scale precisely exponentially P-a.s. Replacing / by the function 5 = 1, it is still 
true that the growth is not precisely exponential - this is readily seen in Claim 
[2] and its proof. 

Since the process in expectation is related to the semigroup S, therefore 
purely exponential scaling indicates the same type of scaling for the expectation 
of the process (locally or globally) . Although, if one is interested in the scaling 
of the process itself (not just the expectation), the case when there is an ad- 
ditional subexponential factor is much harder, in our model the randomization 
of the branching rate f3 helps, as /3 has some 'nice' properties for almost every 
environment u, i.e. the 'irregular' branching rates 'sit in the P-zero set'. 

Define the average growth rate by 



rt = rt{uj) 



log\ZtiLo)\ 



Replace now |^t(w)| by its expectation Zt := E'^\Zt{uj)\ and define 

rt = rt{u;) — - — . 

Recall from Lemma [2l that on a set of full P-measure, 
hm (logt)2/'^(ft-/32) = -c(d,z.). 

t — *oo 

It turns out that an analogous statement holds for rt itself. 
Theorem 12 (LLN). On a set of full P-measure, 

lim {logt)'^^'^{rt - (32) = ~c{d, v) in P'^ - probability. 



(31) 



(32) 



One interprets Theorem [T^ as a kind of quenched law of large numbers. 
Loosely speaking, 

rt~ l32-c{d,v){\ogt)-^/'^ ^rt, t ^ cx), 

on a set of full P-measure. 

While the lower estimate is hard, the upper estimate follows easily from the 
expectation result. We only show here the latter one, but we will explain the 
'philosophy' for the lower estimate too. 

Upper estimate. Let e > 0. Using the Markov inequality along with the 
expectation formula (|?7|) . we have that on a set of full P-measure: 



= P' 

< E 
= exp 



{\ogtf'''{rt~(32)+c{d,y)>e 



{I 



Zt\ > exp 



t - c{d,,y){\ogt)-'/'' + e(logt)-2/'^)] } 



\Zt\ ■ (exp [t (/32 - c{d,iy){logt)-^/^ + eilogt)-^/^ 
et(logi)-2/rf + o(t(logt)-2/rfj 



as t 



□ 
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Lower estimate (sketch). Let us just outline tlie strategy of the proof. 

A key step is introducing three different time scales, £{t), m{t) and t where 
£{t) = o{m{t)) and m{t) = o{t) as t ^ oo. For the first, shortest time interval, 
one uses that there are 'many' particles produced and they are not moving 'too 
far away', for the second (of length m(i) — £{t)) one uses that one particle moves 
into a clearing of a certain size at a certain distance, and in the third one (of 
length t — m{t)) one uses that there is a branching tree emanating from that 
particle so that a certain proportion of particles of that tree stay in the clearing 
with probability tending to one. 

It turns out that for example the following choices of £ and m are appropriate: 
let i{t) and m{t) be arbitrarily defined for t G [0,e], and 

e{t) := ti-V(iogiogt)^ ._ ^1-1/(2 log log t)^ fQj. ^ > > e. 

Also, regarding the Poissonian environment one needs the following fact: let 

]j c{d,u) \uujdj 

(where is the principal Dirichlet eigenvalue of — 5A on the d-dimensional 
unit ball, and ojd is the volume of that ball) and let 

p{e) := M^ogtf/'^ - (loglog^)2, t > 0. 

Then, 

P(3 io{uj) > such that Ml > ^o(w) 3 clearing B{xo,p{i)) with |xo| <l) = l. 
The full proof is beyond the scope of these notes. 



9.3. Some results on the spatial spread 



One may wonder how much the speed (spatial spread) of free BBM reduces due 
to the presence of the mild obstacle configuration. Note that wc are not talking 
about the bulk of the population (or the 'shape') but rather about individual 
particles traveling to very large distances from the origin. 

As is well known, ordinary 'free' branching Brownian motion with constant 
branching rate /32 > has radial speed Let Nt denote the population size 

at t > and let ^fc (1 < A; < Nt) denote the position of the fcth particle (with 
arbitrary labeling) in the population. Furthermore, let m{t) denote a number 
for which u{t,m{t)) 



i, where 



u{t,x) := P 



\\^k{t)\\ < X 

l<k<Nt 



In a classic paper, Bramson considered the one dimensional case and proved 
that as f — > 00, 

m{t) = t^/2^ ^logt + 0(l). (33) 
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Since the one-dimensional projection of a d-dimensional branching Brownian 
motion is a one-dimensional branching Brownian motion, we can utilize Bram- 
son's result for the higher dimensional cases too. Namely, it is clear, that in 
high dimension the spread is at least as quick as in (j33p . The asymptotics (j33p 
is derived for the case P2 — the general result can be obtained similarly. 
Studying the function u has significance in analysis too as u solves 

-wxx + /32(u -u), (34) 



dt 2 



with initial condition 



w(0,-) = l[o,oo)(-)- (35) 

We will see below that the branching Brownian motion with mild obstacles 
spreads less quickly than ordinary branching Brownian motion by giving an 
upper estimate on its speed. 

Remark 13. A related result was obtained earlier by Lee-Torcaso [10], but, 
unlike in (I33p . only up to the linear term and moreover, for random walks 
instead of Brownian motions. Their approach was to consider the problem as 
the description of wave-front propagation for a random KPP equation. o 

Before turning to the upper estimate, we discuss the lower estimate. We are 
going to show that, if in our model Brownian motion is replaced by Brownian 
motion with constant drift 7 in a given direction, then any fixed nonempty ball 
is recharged infinitely often with positive probability, as long as the drift satisfies 

I7I < VW^- _ 

For simplicity, assume that d = I (the general case is similar). Fix the envi- 
ronment (jj. Recall Doob's /i-transform of second order elliptic operators: 

L^{-) ■■=\L[h-). 

Applying an /i-transform with h(x) :— exp(— 7a;) , a straightforward computation 
shows that the operator 

L - -—+ —+f3 
2 dx^ ^ dx ^ 

transforms into 

Then one can show that the generalized principal eigenvalue for this latter op- 

2 

erator is — -f /32 for almost every environment. Since the generalized prin- 

2 

cipal eigenvalue is invariant under /i-transforms, it follows that — -I- /32 > 
is the generalized principal eigenvalue of L. By subsection 13.21 then, any fixed 
nonempty interval is recharged infinitely often with positive probability. 

Turning back to our original setting, the application of the 'spme '-technology 
seems also promising. As we have seen in section |3l, the spine method uses 
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a Girsanov-type change of measure on the branching diffusion. The question 

of uniform integrability for the martingale density can be checked and under 
the new measure, the process can be represented 'directly' by a 'spine' type 
construction. The 'spine' represents a special ^devianf particle performing a 
motion which has important properties from the point of view of the problem 
regarding the original branching process. The question regarding 'going back to 
the original measure' is then related to the UI property of the martingale. 

In our case perhaps it is not very difficult to show the existence of a 'spine' 
particle (under a martingale-change of measure) that has drift 7 as long as 
|7l < VWi- 

We close this section with an upper estimate (without proof) in which the 
order of the correction term is Zander than the 0{logt) term appearing in Bram- 

son's result, namely it is O (^ji^^^^T/d^ ■ (AH orders are meant for t 00.) The 
following theorem says that, whatever (3i G {0,p2) is, loosely speaking, at time 
t the spread is not more than 

Recall the notation: c{d, v) := (^^)^^'^, is the volume of the unit ball in 
R"^, and is the principal Dirichlet eigenvalue of —^A on it. 

Theorem 13. Define the functions 
Then 

nit) = tVm - cid, .)/f . + O (^^) ; and ^f 

:= {no particle has left the n(t) — ball up to t} 
J y supp(Z,)CB(0,n(t)) L 

[o<s<t J 

then on a set of full P -measure, 

liminfP'^(At) > 0. 

9.4- Branching L -diffusion with mild obstacles 

We now generalize the setting for the case when the underlying motion is a 
diffusion. Let P be as before but replace the Brownian motion by an L-diffusion 
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Y on R'*, where L is a second order elliptic operator as in Subsection 12.11 The 
branching L-difFusion with the Poissonian obstacles can be defined analogously 
to the case of Brownian motion. 

The following result demonstrates that the local behavior of the process ex- 
hibits a dichotomy. The crossover is given in terms of the local branching rate 
(32 and Xc{L): as already explained in Subsection 12.11 local extinction occurs 
when the branching rate inside the 'free region' K'^ is not sufficiently large to 
compensate the transience of the underlying L-diffusion; if it is strong enough, 
then local mass grows exponentially. 

There is, however, an interesting new feature of the result: neither (3i nor the 
intensity v of the obstacles play role . 

Claim 1 (Quenched exponential growth/local extinction). Given the environ- 
ment Lo, denote by the law of the branching L-diffusion. 

(i) Let P2 > ~Ac(i) and let v > Q and j3i G (0, /?2) be arbitrary. Then the 
following holds on a set of full V -measure: For any e > and any bounded 
open set%^ B d , 

P'^ ( limsupe("^=-^<=(^)+')*Zf(B) = 00 ) > 0. 

V (Too / 



limsupe(-*-^'^(^»*Zt(B) < 00 = 1. 

(a) Let (32 < — Ac(i) and let v > Q and f3i G (0,/32) be arbitrary. Then 
the following holds on a set of full P-measure: For any bounded open set 
B C K'' there exists a P'^-a.s. finite random time to — tQ{B) such that 
Xt{B) — for all t > to, (local extinction). 

Proof. In order to be able to use Theorem[31 we compare the rate (3 with another, 
smooth (i.e. C^) function V. Recalling that K = := Ua;i6supp(w) -^(^i' °^)' 
let us enlarge the obstacles: 

K*=K*^:^ y P(x„2a). 

a:igsupp(cj) 

Then {K*Y C K"". Recall that P{x) := /3il/f (a;)+/32lif<=(a;) < (32 and let V e 
(7 e (0,1]) witlil!! 

h^K'Y <V<I3. (36) 

In fact it is easy to see the existence of such functions which are even C°° by 
writing (3 ^ (32 — {P2 — (3i)tK and considering the function V :— P2 — {P2 ~ (3i)f, 
where / > Ik* and / is a C"^-function obtained as follows. Let / be a sum 
of compactly supported C°°-functions /„, n > 1, with disjoint support, where 



-"^^The existence of a continuous function satisfying II36I I would of course immediately follow 
from Uryson's Lemma. 
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supp(/„) is in the e„-neighborhood of the nth connected component of K* ^ 
with appropriately small < e„'s. Consider the operator L + V on M.'^ and let 
Ac = Ac(w) denote its generalized principal eigenvalue. Since V G C^, we are in 
the setting of Theorem [3l and furthermore, since V < (32, Xc < \c{L) + /32 for 
every uj. 

On the other hand, one gets a lower estimate on Ac as follows. Fix i? > 0. 
Since /32l(i<'*)c < V , by the homogeneity of the Poisson point process, for almost 
every environment the set {x G M"* | V{x) = P2} contains a clearing of radius R. 
Hence, by comparison, Ac > A'^''^-' , where A'-''^-' is the principal Dirichlet eigenvalue 
of L + /32 on a ball of radius R. Since R can be chosen arbitrarily large and since 
lim/{|oo A'^' = Xc{L) + (32: we conclude that Ac > Xc{L) + (32 for almost every 
environment. 

From the lower and upper estimates, we obtain that 

Ac = Xc{L) + (32 for a.e. u. (37) 

Consider now the branching processes with the same motion component L but 
with rate V , respectively constant rate ^2- The statements (i) and (ii) of Claim 
[T]are true for these two processes by (|37|) and Theorem [3l As far as the original 
process (with rate P) is concerned, (i) and (ii) of Claim [T] now follow by w-wise 
comparison. □ 



References 

[1] S. Albeverio and L.V. Bogachev. (2000) Branching random walk in a cat- 
alytic medium. I. Basic equations. Positivity, 4, 41-100. [MR1740207I 

[2] M. D. Bramson (1978) Minimal displacement of branching random walk. 
Z. Wahrsch. Verw. Gebiete 45(2), 89-108. MR0510529 

[3] M. D. Bramson (1978) Maximal displacement of branching Brownian mo- 
tion. Comm. Pure Appl. Math. 31(5), 531-581. MR0494541 

[4] R. S. Cantrell, C. Cosner (2003) Spatial ecology via reaction-diffusion equa- 
tions. Wiley Series in Mathematical and Computational Biology. John Wi- 
ley and Sons. MR2191264 

[5] Englander, J. (2007) Quenched Law of Large numbers for Branching Brow- 
nian motion in a random medium. Ann. Inst. H. Poincare Probab. Statist., 
to appear 

[6] M. Freidlin, (1985) Functional integration and partial differential equa- 
tions. Annals of Mathematics Studies, 109, Princeton University Press. 
MR0833742 

[7] J. W. Harris, S. C. Harris and A. E. Kyprianou (2006), Further proba- 
bilistic analysis of the Fisher-Kolmogorov-Petrovskii-Piscounov equation: 
one sided travelling- waves, Ann. Inst. H. Poincare Probab. Statist. 42(1), 
125-145. MR2196975 

[8] H. Kesten and V. Sidoravicius (2003) Branching random walk with cata- 
lysts. Electron. J. Probab. 8(5) (electronic). MR 1961167i 



Jdnos Engldnder / Branching diffusions, superdiffusions and random media 



356 



[9] A. E. Kyprianou (2005) Asymptotic radial speed of the support of su- 
percritical branching and super-Brownian motion in R"^. Markov Process. 
Related Fields. 11(1), 145-156. MR2099406 

[10] T. Y. Lee and F. Torcaso (1998) Wave propagation in a lattice KPP equa- 
tion in random media, Ann. Probab., 26(3), 1179-1197. IMR16 34418 

[11] N. M. Shnerb, Y. Louzoun, E. Bettelheim and S. Solomon (2000) The 
importance of being discrete: Life always wins on the surface, Proc. Nat. 
Acad. Sciences 97, 10322-10324. 

[12] N. M. Shnerb, E. Bettelheim, Y. Louzoun, O. Agam and S. Solomon (2001) 
Adaptation of auto catalytic fluctuations to diffusive noise, Phys. Rev. E 63, 
21103-21108. 

[13] J. Xin (2000) Front propagation in heterogeneous media, SIAM Rev. 42(2), 
161-230 (electronic) MR l 778352. 



Jdnos Engldnder/ Branching diffusions, superdiffusions and random media 



357 



10. Generalizations and open problems for mild obstacles 

In this section we suggest some further problems and directions for research. 

10.1. More general branching 

It should also be investigated, what happens when the dyadic branching law is 
replaced by a general one (but the random branching rate is as before). In a 

more sophisticated population model, particles can also die — then the obstacles 
do not necessarily reduce the population size as they sometimes prevent death. 

10.1.1. Supercritical branching 

When the offspring distribution is supercritical, the method of our paper seems 
to work, although when the offspring number can also be zero, one has to con- 
dition on survival for getting the asymptotic behavior. 

10.1.2. (Sub) critical branching 

Critical branching requires an approach very different from the supercritical 
one, since taking expectations now does not provide a clue: E'^\Zt{uj)\ = 1, 
Vi > 0, Vw e O. 

Having the obstacles, the first question is whether it is still true that 

P'^ (extinction) = 1 Vw e O. 

The answer is yes. To see this, note that since |^| is still a martingale, it has 
a nonnegative a.s. limit. This limit must be zero; otherwise \Z\ would stabilize 
at a positive integer. This, however is impossible because following one Brow- 
nian particle it is obvious that this particle experiences branching events for 
arbitrarily large times. 

Setting /3i = 0, the previous argument still goes through for almost all en- 
vironments, because for almost all environments, the generic Brownian particle 
will eventually branch, and so in fact it will branch infinitely often. Let r denote 
the almost surely finite extinction time for this case. One of the basic questions is 
the decay rate for P'^(r > t). Will the tail be significantly heavier than 0{l/i)l 
[Of course 1/t would be the rate without obstacles.] 

The subcritical case can be treated in a similar fashion. In particular, the 
total mass is a supermartingale and P'^ (extinction) = 1 Vw e fi. 

10.2. Superprocesses with mild obstacles 

An alternative view on the BBM with mild obstacles is as follows. Arguably, the 
model can be viewed as a catalytic BBM as well — the catalytic set is then K'^ 
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(in the sense that branching is 'intensified' there). Catalytic spatial branching 
(mostly for superprocesses though) has been the subject of vigorous research 
in the last twenty years initiated by Dawson, Fleischmann and others. In those 
models the individual branching rates of particles moving in space depend on 
the amount of contact between the particle ('reactant') and a certain random 
medium called the catalyst. The random medium is usually assumed to be a 
'thin' random set (that could even be just one point) or another superprocess. 
Sometimes 'mutually' or even 'cyclically' catalytic branching is considered. (See 
the next section for more on catalytic branching.) 

Although there is a large amount of ongoing research on catalytic superpro- 
cesses, in those models, one usually cannot derive sharp quantitative results. 

In light of this connection with catalytic superprocesses, a further goal is to 
generalize our setting by defining superprocesses with mild obstacles analogously 
to the BBM with mild obstacles. 

Consider the (L, /3, a, R'^)-superdiffusion, X. The definition of the superpro- 
cess with mild obstacles is straightforward: the parameter a on the (random) 
set K is smaller than elsewhere. 

Similarly, one can consider the case when instead of a, the 'mass creation 
term' /3 is random, for example with /3 defined in the same way (or with a 
mollified version) as for the discreet branching particle system. Denote now by 
P'^ the law of this latter superprocess for given environment. We suspect that 
the superprocess with mild obstacles behaves similarly to the discrete branch- 
ing process with mild obstacles when Xc{L + (3) > and /*"(•) is replaced by 
P'^(- I X survives). The upper estimate can be carried out in a manner similar 
to the discrete particle system, as the expectation formula is still in force for 
superprocesses. 

10.3. Unique dominant branch 

Having discussed the growth rate of the population, the next step can be as fol- 
lows. Once one knows the global population size \Zt\, the model can be rescaled 
(normalized) by \Zt\, giving a population of fixed weight. In other words, one 
considers the discrete probability measure valued process 

Then the question of the shape of the population for Z for large times is given 
by the limiting behavior of the random probability measures Zt^ t > Q. (Of 
course, not only the particle mass has to be scaled, but also the spatial scales 
are interesting — see last paragraph.) 

Can one for example locate a unique dominant branch (or, even if it is not 
unique, there are perhaps just a 'few' of them) for almost every environment, 
so that the total weight of its complement tends to (as t ^ oo) zero? 

The motivation for this question comes from the proof of the lower estimate 
for Theorem [121 It seems conceivable that for large times the 'bulk' of the 
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population will live in a clearing(s) within distance £{t) and with radius 
p{t) :-i?o[log^W]i/'^-[loglog^(t)]2. 

10.4- Strong Law 

Last but not least, one cannot be fully satisfied with the convergence in proba- 
bility in Theorem[3l To prove a.s. convergence (which, if true, is a kind of SLLN) 
seems to be rather challenging. In particular, the easy upper estimate does not 
work any longer. 
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11. Some catalytic branching problems 

11.1. A discrete model 

Let us start with a discrete catalytic model studied by H. Kesten and V. Sido- 
ravicius. Here the branching particle system on X'^ is so that its branching is 
catalyzed by another autonomous particle system on lA. There are two types 
of particles, the A-particles ('catalyst') and the -B-particles ('reactant'). They 
move, branch and interact in the following way. Let Na{x, s) and Nb{x, s) de- 
note the number of A- [resp. i3-]particles at a; e Z'' and at time s € [0, oo). (All 
Na{x, 0) and Nb{x, 0), x € Z"^ are independent Poisson variables with mean fiA 
(mb).) 

Every A-particle (i?-particle) performs independently a continuous-time ran- 
dom walk with jump rate Da {Db)- In addition a -B- particle dies at rate r, 
and, when present at x at time s, it splits into two in the next ds time with 
probability PNa{x, s)ds + o{ds). Conditionally on the system of the A-particles, 
the jumps, deaths and splitting of the B-particles are independent. 

Kesten and Sidoravicius showed for instance that for large P there exists a 
critical Tc such that it separates local extinction regime (r > Vc) from local 
survival regime (r < Tc). They also explained why the B-particles may survive 
if P is large enough compared to the other parameters, even if one starts with 
only one A-particle and one B-particle in the entire system. 

11.2. Catalytic superprocesses 

The definition of catalytic superprocesses requires more care. The simplest case 
was introduced by D. Dawson and K. Fleischmann who investigated a one 
dimensional critical super-Brownian motion with a single point catalyst. This 
means formally that we set [3 — Q and a = 5q (the Dirac-delta at zero). In- 
tuitively this means that branc hing occurs only at a single point but there it 
occurs at an 'infinitely high' ratqlj. The precise definition is given through the 
log-Laplace equation, in which case ([4]) is replaced by an appropriate integral 
equation. 

One of the interesting properties of this process is that if P denotes the 
corresponding probability, then 

P^(|Xt| > 0, > 0, but \Xt\ as t oo) = 1 

for aU ^ /i G A^/, and Xt{B) ^ in -probability for aU B CC M, even 
if /i is the Lebesgue measure. (The definition of a superprocess starting with a 
u-finite measure is straightforward: it is the independent sum of superprocesses 
starting with finite measures.) 

Another interesting case is, when, instead of the singular coefficient 5q in 
the branching mechanism, one imagines that the catalytic set, where branching 



Recall that since /3 = 0, one can now identify a with the 'clock' 
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occurs, is given by the support of another superprocess (in low dimensions). Note 
that this is aheady the case of a random catalyst, that is, we are talking about 
a superprocess in a random medium (but, unlike in the cases we discussed, the 
medium evolves in time.) Since branching is critical, it is detrimental for the 
reactant to spend too much time in catalytic areas. 

For d < 3, Dawson and Fleischmann investigated the qualitative properties 
of this superprocess. They showed for example, that in dimension one, if both 
the catalyst and the reactant starts from Lebesgue measure then the reactant 
process survives locallyEl This is dramatically different behavior from that of a 
classical one dimensional super-Brownian motion, since the latter would suffer 
local extinction even starting from fi = Lebesgue. 

In two dimension, again starting both processes from Lebesgue measure, Daw- 
son and Fleischmann proved that although the catalyst dies out locally, this is 
true only in probability. In particular, there is not a finite time after which a 
given ball becomes and remains empty. Consequently, at late times T, huge clus- 
ters of the catalyst (their height turns out to be of order log T ) come back to 
any given finite window in and (since critical binary branching with infinite 
rate degenerates to pure killing) could kill all the (recurrent) reactant particles 
there. Thus it is somewhat surprising that in fact, the reactant survives locally. 

Now let us go one step further, and imagine that we have two one dimensional, 
critical super-Brownian motions and they catalyze each other. That is, the first 
process only branches in the presence of the second process, and vice versa 
{'mutually catalytic branching') . Then, of course it becomes a nontrivial business 
to rigorously prove that such a model is well defined; in fact the log-Laplace 
equation approach breaks down. Such a model was first studied by D. Dawson 
and E. Perkins [7]. The definition uses a system of stochastic partial differential 
equations (SPDE's) for the densities of the two superprocesses, which by a result 
of L. Mytnik [19], uniquely determine the model. 

Further models, including the cyclically catalytic setting, are beyond the 
scope of these lecture notes. We suffice here by merely defining the latter model. 
Cyclically catalytic super-Brownian motion is a spatial branching process with 
K types of populations involved {K >2). Each process X'' is (catalytic) super- 
Brownian motion in one dimension but with a branching rate that is governed by 
the density of the succeeding type X^^^. Hence the K single processes interact 
in a cyclic fashion: 

0->K-l->K-2-> >1^0. 

The literature of catalytic spatial branching processes and superprocesses is 
huge and we have just glimpsed at it. The interested reader is referred to the 
survey paper [16] . 



^^In fact, they showed that for almost all realizations of the catalytic medium (starting 
with Lebesgue measure), the reactant at t converges (stochastically) to the starting Lebesgue 
measure. 
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12. Concluding remarks 

Spatial branching processes is a huge area. So is spatial random media. Although 
we tried to cull some interesting problems, there is no doubt that many papers 
that should have been mentioned wore unfortunately skipped. 

As far as superprocesses are concerned, it is impossible even to sketch the 
research that has been carried out by Aldous, Dawson, Delmas, Donelly, Dynkin, 
Evans, Flcischmann, Klcnkc, Kurtz, Lc Gall, Lc Jan, Mortcrs, Mytnik, Mueller, 
Perkins, Pinsky and others. (And for discrete branching processes, there are 
those work of Athreya, Biggins, S. Harris, T. Harris, Jagers, Kyprianou, Ney, 
Pcmantlc, Wakolbingcr, and the whole French school, to name just a few...) 

Similarly, we did not even attempt to describe all the contributions made to 
the field of spatial random media by scholars like Bolthausen, Comets, Cranston, 
Donskcr, Gartner, den Hollander, van den Berg, Molchanov, Sznitman, Varad- 
han, Zeitouni, Zerner and others. 

My hope is that those whose work was not mentioned in these notes will 
(forgive the author and) send me their critical comments. 



